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16.  Abifscl 

v  New  computational  methods  have  been  applied  in  an  investigation  of  the 
nonlinear  surge  motions  and  springing  loads  in  the  tendons  of  tension  leg 
platforms  (TLPs).  The  second-order  wave  drift  forces  as  well  as  the 
linear  hydrodynamic  coefficienta  and  the  wave-induced  exciting  forces  are 
computed  by  a  three-dimensional  hybrid-finite-element  method  (HFEM). 

Both  a  new  formula  for  predicting  the  viscous  drag  forces  and  the 
conventional  Morison  drag  formula  have  been  used  in  a  time-domain 
computational  procedure  .to  predict  the  nonlinear  surge  motions  for  two 
TLP  configurations  due  to  wind,  current  and  wave  excitations.  The 
wave-induces  first-and  second-order  pitch  springing  exciting  moments  and 
the  resulting  springing  tension . loads  in  the  tendons  are  computed  for 
regular  waves,  two  waves  (wave  groups)  and  irregular  waves  using  a 
short-wave  approximation  method.  Upper  bound  estimates  of  the  surge 
siotions  and  the  tendon  springing  loads  are  computed  for  some  selected 
extreme  environmental  conditions.  "Average"  surge  displacements  and 
"average"  springing  loads  are  also  computed  for  some  selected  typical 
operational  conditions.  Finally,  the  surge  motions  and  the  springing 
loads  of  a  TLP  with  one  tendon  removed  have  been  predicted  and  compared 
to  the  results  for  the  same  TLP  with' all  of  the  tendons  intact. 
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ABSTRACT 


New  computational  qjefhodfi  have  been  applied  ip  an  investigation  of 
the  nonlinear  surge  motions  and  springing  loads  in  the  tendons  of  tension 
leg  platforms  (TIPs).  The  second-order  wave  drift  forces  as  well  as  the 
linear  hydrodynamic  coefficients  and  the  wave-induced  exciting  forces  are 
computed  by  a  three-dimensional  hybrid-finite-element  method  (HFEM).  Both 
a  new  formula  for  predicting  the  viscous  drag  forces  and  the  conventional 
Morison  drag  formula  have  been  used  in  a  time-domain  computational  procedure 
to  predict  the  nonlinear  surge  motions  for  two  TLP  configurations  due  to 
wind,  current  and  wave  excitations.  The  wave-induced  first-  and  second-order 
pitch  springing  exciting  moments  and  the  resulting  springing  tension  loads  in 
the  tendons  are  computed  for  regular  waves,  two  waves  (wave  groups)  and  ir¬ 
regular  waves  using  a  short-wave  approximation  method.  Upper  bound  estimates 
of  the  surge  motions  and  the  tendon  springing  loads  are  computed  for  some 
selected  extreme  environmental  conditions.  "Average"  surge  displacements 
and  "average"  springing  loads  are  also  computed  for  some  selected  typical 
operational  conditions.  Finally,  the  surge  motions  and  the  springing  loads 
of  a  TLP  with  one  tendon  removed  have  been  predicted  and  compared  to  the 
results  for  the  same  TLP  with  all  of  the  tendons  intact. 
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INTRODUCTION 


1 . 1  Objectives 

The  main  objectives  of  this  project  have  been  to  investigate  the 
following  nonlinear  TLP  responses: 

•  surge  motions  due  to  nonlinear  wave-induced  drift  forces,  wind 
loading  and  viscous  interactions  between  current  and  oscillatory 
motions 

§  additional  tension  loads  in  the  tendons  due  to  nonlinear  wave- 
induced  springing  loads,  and 

•  motions  and  loads  for  reduced  number  of  tendons. 

1.2  Background 

1.2.1  Tension  Leg  Platform  Concept 

In  order  to  minimize  the  wave-induced  tension  loads  and  the  wave- 
induced  surge  displacements,  TLPs  are  designed  in  such  a  way  that  the  natural 
periods  of  oscillations  for  the  vertical  heave  and  pitch  motions  and  for 
the  horizontal  surge  and  yaw  motions  are  far  outside  the  range  in  which  the 
wave  energy  of  the  ocean  is  a  maximum.  TLPs  usually  have  a  natural  period 
in  heave  of  about  2-3  seconds  and  a  natural  period  in  surge  of  about  100- 
140  seconds. 

When  waves  pass  a  floating  body,  it  will  be  excited  by  two  types 

of  forces: 

•  The  primary  linear  exciting  forces  which  have  the  same  period  as 
the  wave  period,  and 
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•  The  secondary  nonlinear  exciting  forces  which  have  steady  and 
oscillatory  components  with  periods  far  outside  the  range  of 
the  typical  wave  periods.  These  nonlinear  forces  are  often 
referred  to  as  the  slowly-varying  drift  and  high-frequency  spring¬ 
ing  exciting  forces. 

An  important  design  aspect  of  the  TLP  is  that  the  nonlinear  springing  and 
drift  forces  may  introduce  excitation  at  the  natural  heave  and  surge  fre¬ 
quencies  of  the  TLP.  Therefore,  it  is  possible  that  the  high-frequency 
springing  forces  can  cause  additional  tension  loads  in  the  tendons  and  that 
the  slowly- varying  drift  forces  can  cause  surge  displacements  which  are  as 
large  if  not  larger  than  the  loads  and  displacements  caused  by  the  primary 
linear  wave- induced  exciting  forces. 

1.2.2  Prediction  of  TLP  Responses 

It  is  relatively  easy  to  predict  the  primary  linear  exciting 
forces  and  the  resulting  motions  and  tension  loads.  Since  these  primary 
motions  occur  at  frequencies  far  removed  from  the  TLPs  natural  frequencies, 
they  can  be  predicted  accurately  using  a  frequency  domain  approach  with 
linear  coupled  differential  equations  ignoring  the  nonlinear  viscous  damp¬ 
ing. 

On  the  other  hand,  the  secondary  motions  and  loads  caused  by  the 
nonlinear  wave-induced  excitation  is  very  difficult  to  predict  accurately. 
The  major  problem  areas  are  the  three  following: 

•  Accurate  prediction  of  the  nonlinear  drift  and  springing  excit¬ 
ing  forces  in  regular  waves.  (Prediction  methods  only  exist  for 
the  second-order  forces . ) 

•  Accurate  prediction  of  the  slowly-varying  surge  responses  and  the 
high-frequency  springing  responses  which  both  occur  at  the  natural 
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TLP  frequencies  and  hence  are  governed  by  the  nonlinear  viscous 
damping  forces.  (The  accuracy  of  existing  methods  for  predict¬ 
ing  the  viscous  damping  for  TLPs  is  questionable.) 

•  Accurate  predictions  of  drift  and  springing  responses  in  irregu¬ 
lar  seas  is  a  very  complicated  problem.  Not  only  is  the  viscous 
damping  nonlinear  with  complicated  interaction  effects  between 
the  steady  current  and  the  oscillatory  motions,  but  more  impor¬ 
tantly,  the  nonlinear  forces  in  irregular  seas  is  due  to  non¬ 
linear  interactions  between  the  different  wave  components. 

1.2.3  The  TLP  Motion  and  Load  Computer  Code  System 

A  large-scale  computer  code  system  has  been  developed  by  SAI  specif¬ 
ically  for  the  purpose  of  analyzing  the  nonlinear  dynamics  of  TLPs.  Figure 
1-1  shows  a  schematic  layout  of  the  TLP  Motion  and  Load  Computer  Code  Sys¬ 
tem.  The  system  consists  of  a  total  of  seven  individual  computer  codes 
which  are  grouped  in  three  major  parts: 

Part  I:  Potential-Flow  Predictions  by  a  Hybrid-Finite-Element 
Method  (HFEM)  Computer  Program 

Part  II:  Surge-Response  Predictions  by  Time-Domain  Integrations 
Method 

Part  III:  Springing-Response  Predictions  by  a  New  High-Frequency 
Method 

This  computer  code  system  is  a  result  of  several  man-years  of  research 
work  and  code  development.  It  is  probably  one  of  the  most  advanced  TLP 
dynamic  response  prediction  systems  presently  in  existence. 
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Figure  1-1.  TLP  '-lotions  and  Loads  Computer  Code  System 


TENDON  LOADS 
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Summary  of  Resul ts 

Surge  motions  and  tendon  tension  springing  loads  have  been  com¬ 
puted  for  two  TLPs.  TLP  HI  is  a  2000-foot  depth  design  and  TLP  HZ  is  a 
3000-foot  depth  design.  More  detailed  descriptions  of  their  geometry  are 
given  in  Appendix  A.  The  main  underwater  structure  of  both  of  the  platforms 
consists  of  four  vertical  corner  legs  and  four  horizontal  pontoons.  Both 
platforms  have  four  vertical  tendons  at  each  of  the  four  legs. 

Motions  and  loads  have  been  predicted  for  extreme  environmental 
conditions  as  well  as  for  a  typical  operational  condition. 

1.3.1  Surge  Motion  Results 

The  surge  motions  were  computed  by  the  nonlinear  time-domain  surge 
motion  code  for  a  large  number  of  wave,  wind  and  current  conditions.  The 
most  important  findings  are  as  follows: 

(a)  Upper  bound  estimates 

For  TLP  HI  it  was  found  that  for  the  environmental  conditions 
considered  here  the  maximum  surge  displacement  is  152.3  feet  (166.6 
feet)*.  This  maximum  surge  displacement  will  occur  in  a  condition 
consisting  of  maximum  unsteady  wind,  maximum  steady  current  and 
regular  waves  with  period  T  =  9.5  seconds.  This  is  the  wave  period 
for  which  the  steady  wave-induced  drift  force  is  maximum.  It  is 
here  assumed  that  the  wave  steepness  measured  as  the  ratio  between 
the  wave  height  and  the  wave  length,  H/A,  is  0.10. 


*  In  this  section,  the  results  given  in  parentheses  have  been  obtained  by 
the  conventional  Morison  drag  formula,  whereas  the  results  not  in  paren¬ 
theses  have  been  obtained  by  a  new  drag  formula. 


1-5 


It  should  be  emphasized  that  the  surge  displacement  is  smaller 
for  the  maximum  design  wave  condition  (H  -  80  feet  and  T  -  14.0 
seconds)  than  it  is  for  the  condition  stated  above.  Our  computa¬ 
tions  show  that  the  maximum  surge  displacement  is  102.4  feet 
(127.4  feet)  for  the  maximum  design  wave  condition  with  maximum 
current  and  maximum  unsteady  wind. 

These  are  extreme  upper  bound  estimates  which  will  only  occur 
under  extreme  and  rare  conditions.  Further  investigations  of  the 
probability  of  the  occurrance  of  such  rare  events  is  needed. 

(b)  Typical  operational  estimates 

The  surge  motions  have  also  been  computed  for  a  "typical" 
irregular  wave  condition  with  significant  wave  height, 
equal  to  10  feet  and  with  steady  current  and  steady  wind  of  a 
strength  typical  for  this  sea  condition.  For  TLP  # 1  the  total 
"average"  surge  displacement  is  estimated  to  be  about  26  feet. 

Note  that  the  estimates  show  that  the  amplitudes  of  the  wave 
frequency  and  the  slowly-varying  components  are  less  than  two  and 
four  feet  respectively  and  that  the  steady  surge  can  be  expected 
to  be  about  20  feet. 

Also  it  should  be  pointed  out  that  the  investigation 
revealed  a  large  difference  between  the  surge  results  obtained 
using  the  two  viscous  drag  formulas.  The  large  difference  between 
the  two  formulas  is  due  to  the  fact  that  the  conventional  Mori  son 
quadratic  formula  predicts  a  very  large  additional  steady  displace¬ 
ment  due  to  interaction  effects  between  the  wave  frequency  compo¬ 
nents  and  the  steady  current.  It  is  believed  that  this  interac¬ 
tion  effect  is  overpredicted  by  the  conventional  quadratic  drag 
formula. 


such  large  slowly-varying  motions,  we  need  to  know  more  about  the  probabil¬ 
ity  and  statistics  of  such  events  (or  events  close  to  these  events),  for 
typical  general  irregular  sea  conditions  it  can  be  expected  that  the  slowly- 
varying  motions  are  quite  small,  probably  less  than  say  4  feet. 


The  steady  surge  displacements,  on  the  other  hand,  can  be  expected 
to  be  much  larger  than  the  unsteady  motions.  For  some  of  the  conditions 
investigated,  we  have  the  following  results  for  the  steady  surge. 

Steady  Surge 


•  Maximum  current 

(a)  with  zero  unsteady  surge  14.4 

(b)  with  maximum  wave-freuqency  surge  42.3* 


•  Maximum  steady  wind 


52.5 


•  Steady  wave- induced  drift 

(a)  regular  wave  with  maximum  drift  64.1 

(b)  typical  irregular  sea  -  10')  2.3 


Hence,  for  typical  environmental  conditions  one  may  expect  up  to  about  20 
feet  of  steady  surge.  However,  the  extreme  steady  surge  displacement  may  be 
considerably  larger. 


2.2  HFEM  Potential  Flow  Calculations 


The  three-dimensional  linear  diffraction  analyses  were  performed 
using  the  SAI  Hydro-Finite-Element  Method  (HFEM)  code  utilizing  the  vertical 
(perpendicular)  planes  of  symmetry.  The  HFEM  diffraction  code  is  a  powerful 
numerical  method  for  calculating  the  hydrodynamic  forces  and  the  wave- 
induced  responses  and  loads  of  three-dimensional  offshore  structures.  The 

*  As  predicted  by  Morison  quadratic  drag  formula. 
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These  results  show  that  large  wave-frequency  surge  motions  can  only  be 
expected  for  very  extreme  wave  conditions  as,  for  example,  the  maximum 
design  regular  wave.  For  typical  wave  conditions  it  can  be  expected  that 
the  amplitude  of  the  wave- frequency  component  will  be  less  than  two  feet. 

surge  motions  are  given  below 

Surge  Amplitude 

40.9' 

=  20' 

•  Two  regular  waves  (wave  grouping) 

(a )  zero  current  ~  58 ' 

(b)  maximum  current  ~25‘ 

•  Typical  irregular  waves  (H^  =  10')  2.6' 

These  results  show  that  large  slowly-varying  natural-frequency  surge  motions 
can  be  expected  under  "ideal"  severe  conditions  when  the  excitation  fre¬ 
quency  is  precisely  equal  to  the  natural  frequency.  For  the  unsteady  wind 
case  it  is  assumed  that  the  unsteady  part  of  the  wind  has  one  single  fre¬ 
quency  component  with  amplitude  equal  to  107.  of  the  maximum  steady  wind 
speed  and  frequency  equal  to  the  natural  frequency.  For  the  two-wave  case, 
it  is  assumed  that  the  difference  frequency  for  the  two  wave  components  is 
equal  to  the  natural  frequency: 

oo  2  w  ^  -  ojp  • 

This  implies  that  there  will  be  a  second-order  si owly-varyi ng  wave-exciting 
force  with  excitation  frequency  equal  to  the  difference  frequency.  Both 
of  these  natural-frequency  excitation  conditions  are  very  "idealistic" 
conditions.  Before  we  can  make  any  predictions  about  the  probability  of 


to  these  viscous  interaction  effects  can  typically  reduce  the  amplitude  of 
the  slowly-varying  motions  by  as  much  as  a  factor  of  four. 


The  steady  surge  displacement  is  due  to 

•  steady  wind  forces, 

•  steady  nonlinear  wave- induced  drift  forces, 

•  steady  current  forces,  and 

•  nonlinear  interaction  between  current  and  wave- frequency  motions 

and  slowly-varying  motions. 

In  this  project  we  have  computed  the  surge  displacement  due  to  all  of  the 
above-stated  steady  forces;  however,  of  main  interest  in  this  study  is  the 
computation  of  the  steady  surge  due  to  the  nonlinear  wave-drift  force  and 
the  nonlinear  viscous  interaction. 


We  shall  now  present  some  results  for  the  different  surge  dis¬ 
placement  components  so  that  we  can  get  a  better  feeling  for  the  order  of 
magnitude  of  the  different  components  and  their  relative  importance.  First 
we  shall  consider  the  wave-frequency  surge  displacement  components.  The 
amplitude  of  the  wave-frequency  component  for  three  wave  conditions  is 
given  below. 


•  Maximum  regular  wave 

(H  =  80’,  X  =  1000'  and  T  -  14  sec) 


Surge  Amplitude 

22. r 


§  Steep  regular  wave  with  typical  wave  length 

(H  =  46’ ,  X  =  462'  and  T  =  9.5  sec)  1.7' 

•  Typical  irregular  waves  (H..-  =  10’)  0.9' 


We  shall  now  take  a  look  at  the  different  components  of  the  surge 
displacement.  The  surge  mo Mop  of  a  TLP  can  he  separate^  into  three  main 
parts : 

x(t)  =  X  +  Hx.  cosm.t  +  EX.  cosm.t 
o  j  i  1  j  J  J 


The  first  term,  XQ,  is  the  mean  s teady  displacement.  The  second  term, 

E  X-  cosm-t,  consists  of  first-order  frequency  components  which  have  fre¬ 
quencies  identical  to  the  wave  frequencies.  We  shall  refer  to  these  as  the 
linear  wave-frequency  components.  The  third  and  final  term,  E  x  j  cosiOjt, 
consists  of  frequency  components  which  have  frequencies  near  or  equal  to 
the  natural  surge  frequency.  Note  that  the  natural  frequency  components 
are  far  removed  from  the  wave-frequency  components.  We  shall  refer  to  these 
as  the  nonlinear  slowly- varying  components. 

It  is  important  to  recognize  that  the  linear  wave- frequency  com- 
ponents  are  caused  by  linear  wave-induced  excitation  and  that  the  motions 
at  these  frequencies  are  governed  by  the  inertia  forces  and,  furthermore, 
that  viscous  damping  has  practically  no  effect  or.  these  frequency  components. 
Hence,  the  wave  frequency  motion  can  easily  be  predicted  by  thfe  well  estab¬ 
lished  linear  superposition  approach. 

On  the  other  hand,  the  slowly-varying  components  at  or  near  the 
natural  frequency  are  caused  by  nonlinear  wave-induced  excitation  and/or 
wind  excitation.*  Furthermore,  since  these  motions  are  at  or  near  the  natu¬ 
ral  frequency,  they  are  to  a  large  extent  governed  by  nonlinear  viscous 
damping  forces.  Note  that  the  damping  of  the  slowly-varyi ng  motions  is  very 
much  affected  by  nonlinear  viscous  interactions  between  the  slowly-varying 
motion,  the  steady  current  and  the  wave-frequency  motions.  The  damping  due 

*  Note  that  nonlinear  viscous  interactions  between  the  steady  current  and 
the  wave-frequency  components  may  also  result  in  slowly-varyi ng  excita¬ 
tion  which  has  not  been  considered  in  this  investigation. 


The  results  presented  in  Table  2-1  show  that  for  the  maximum  regu¬ 
lar  wave  condition  with  maximum  steady  wind  and  current,  the  maximum  surge 
displacement  computed  by  the  SAI  drag  formula  is  92.7  feet  and  computed  by 
the  quadratic  drag  formula  is  118.8  feet.  The  large  difference  between  the 
two  drag  formulas  is  due  to  the  fact  that  the  quadratic  formula  predicts  a 
very  large  additional  steady  displacement  due  to  interaction  effects  between 
wave  frequency  and  current.  It  is  believed  that  this  interaction  effect  is 
overpredicted  by  the  quadratic  drag  formula.  Furthermore,  for  the  maximum 
regular  wave  condition  it  is  seen  in  Table  2-1  that  the  unsteady  wind 
increases  the  surge  by  9.7  feet  when  the  SAI  drag  formula  is  used,  whereas 
it  increases  by  only  8.9  feet  when  the  quadratic  drag  formula  is  used. 

The  second  wave  condition  presented  in  Table  2-1  is  the  regular 
wave  condition  which  will  result  in  the  maximum  drift  force.  The  results 
presented  in  Section  2.2  show  that  the  maximum  drift  force  will  occur  for 
a  regular  wave  with  T  =  9.5  seconds.  Our  predictions  show  that  the  surge 
displacement  for  this  wave  condition  with  maximum  steady  wind  and  current 
will  be  133.5  feet  when  using  the  SAI  drag  formula  and  153.8  when  using  the 
quadratic  formula.  The  unsteady  wind  results  show  152.3-foot  surge  displace¬ 
ment  for  the  SAI  drag  case  and  166.6  feet  for  the  quadratic  drag  case.  Note 
that  the  surge  displacements  are  substantially  larger  for  the  maximum  drift 
wave  case  than  for  the  maximum  wave  condition. 

Significant  surge  displacements  are  also  shown  in  Table  2-1  for 
the  typical  irregular  case  as  computed  by  the  SAI  drag  formula.  With  zero 
current  the  surge  displacement  is  only  11.6  feet  and  with  maximum  current 
it  is  25  feet.  Therefore,  for  typical  general  irregular  wave  conditions, 
the  surge  displacements  can  be  expected  to  be  very  small. 

Total  surge  displacements  were  not  computed  for  the  two  regular 
wave  case  (wave  grouping).  For  this  case  only  the  slowly-varying  wave  exci¬ 
tation  and  response  were  investigated  since  these  are  the  most  important 
aspects  of  wave  grouping  excitation. 


•  Typical  irregular  wave  condition  (Hj^  -  10  ft) 

Different  combinations  of  wind  and  current  conditions  have  been  used  with 
these  wave  conditions. 


A  summary  of  the  total  surge  displacements  for  the  maximum  regular 
wave,  the  regular  wave  with  maximum  drift  and  the  typical  irregular  wave 
conditions  is  presented  in  Table  2-1  for  the  TLP  #1  configuration  (see 
Appendix  A  for  TLP  dimensions).  Results  are  shown  as  obtained  by  the  new 
SA1  drag  formula  and  by  the  quadratic  Morison  drag  formula. 

Table  2-1.  Total  surge  displacements  for  TLP  ill  for  two  regular- 

wave  conditions  and  a  typical  irregular  wave  condition. 


Maximum  Surge  Displacement 
( i n  feet) 

SA1  Drag 

Quadratic 

Drag 

!  Maximum  Regular  Wave  Conditions 

Wave:  Maximum  Regular  Wave 
(H  -  80  and  T  =  14.0) 

Current:  Maximum  Steady 

Wind:  (a)  Maximum  Steady 

92.7 

118.8 

(b)  Unsteady 

102.4 

127.7 

Regular  Wave  with  Maximum  Wave  Drift 

Wave:  Regular  Maximum  Drift 
(H  -  46.4  and  T  -  9.5) 

Current:  Maximum  Steady 

Wind:  (a)  Maximum  Steady 

133.5 

153.8 

(b)  Unsteady 

152.3 

166.6 

Typical  Irregular  Wave  Conditions 

Wave:  Irregular  Waves  (H,,n  =  10) 

Wind:  Steady  (39  ft/ sec)  ' 

Current:  (a)  Zero  Current 

11.6* 

(b)  Maximum  Steady 

25.0* 

*  This  is  the  significant  surge  displacement,  x^- 
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the  linear  wave-exciting  forces  and  the  second-order  wave-drift  forces  are 
computed  by  the  hybrid-f ini te-element  method.  These  quantities  are  then 
used  as  input  for  the  nonlinear  time-domain  surge-motion  computations. 

Computations  have  shown  that  three-dimensional  hydrodynamic  inter¬ 
action  effects  between  the  legs  and  the  pontoons  have  an  important  influence 
on-the  values  of  the  predicted  added  mass,  wave  damping  and  exciting  forces. 
Hence,  the  complete  three-dimensional  computation  has  to  be  performed.  Com¬ 
putations  of  hydrodynamic  quantities  that  are  made  for  the  legs  and  for  the 
pontoons  separately  and  then  simply  added  will  not  be  sufficiently  accurate 
to  correctly  predict  the  surge  motion  of  the  TLP. 

The  viscous  forces  are  predicted  by  a  new  drag  formula  as  well  as 
by  the  conventional  quadratic  Morison  drag  formula.  This  new  formula  has 
been  developed  because  it  is  believed  that  the  Morison  drag  formula  over¬ 
predicts  the  magnitude  of  the  viscous  interaction  effects  between  current, 
wave-frequency  motions  and  slowly-varying  motions.  The  surge  motions  due 
to  simultaneous  action  of  the  wind,  current  and  wave  excitation  are  pre¬ 
dicted  by  a  numerical  time-integration  procedure  which  includes  all  of  the 
important  nonlinear  effects. 

2.1.2  Summary  of  Surge  Motion  Results 

In  this  study,  we  have  investigated  the  surge  motions  for  four 
different  wave  conditions: 

•  Maximum  regular-wave  condition  (H  =  80  ft,  X  =  1000  ft  and  T  =  14  sec 

•  Regular  wave  with  maximum  steady  drift  force  (H=46.4  ft,  X  =  462 

ft  and  T  =  9. 5  sec) 


Two  regular  waves  (Wave  grouping)  (T^  =  9.7  sec  and  T^  -  8.9  sec) 
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SURGE  MOTION  INVESTIGATION 


2.1  Introduction 

2.1.1  Objective  and  Approach 

The  main  objective  of  this  part  of  the  project  has  been  to  inves¬ 
tigate  the  surge  motions  due  to  nonlinear  wave-induced  drift  forces,  non¬ 
linear  wind  loading  and  nonlinear  viscous  interactions  between  current  and 
oscillatory  motions.  Two  TLP  conf igurations  are  analyzed  for  selected 
severe  and  typical  average  wind,  wave  and  current  conditions.  An  extensive 
evaluation  and  comparison  of  the  surge  displacements  due  to  the  important 
nonlinear  viscous  interaction  effects  between  current,  wave- frequency  motions 
and  slowly-varying  motions  as  predicted  by  a  new  viscous  force  formula  and 
by  the  conventional  Morison  force  formula  are  presented.  Unfortunately, 
there  are  practically  no  experimental  ersults  available  for  these  important 
nonlinear  viscous  interaction  effects.  There  is  also  a  lack  of  experimental 
data  for  the  nonlinear  surge  motions  under  controlled  environmental  condi¬ 
tions  for  which  computations  can  be  performed.  (See  Appendix  D  for  a  dis¬ 
cussion  of  comparisions  between  theory  and  experiments.) 

A  new  computational  method  developed  by  SAI  (Salvesen,  et  al . , 

1982)  is  here  used  to  investigate  the  nonlinear  surge  motions  for  TLPs.  The 
computations  consist  of  three  main  parts: 

•  potential -flow  computations; 

•  viscous-force  predictions;  and 

•  nonlinear  time-domain  surge  computations. 

The  potential-flow  calculations  are  performed  first  by  a  three-dimensional 
hybrid-finite-element  method.  The  added  mass,  the  wave-damping  coefficients, 


In  this  investigation  we  have  only  considered  additional  natural 
frequency  tendon  loads  due  to  periodic  springing  excitation.  Wave  impact 
excitation  may  also  result  in  large  natural  frequency  tendon  lo.nls  often 
referred  to  as  "ringing"  loads.  It  is  recommended  that  an  investigation  of 
the  wave-impact  "ringing"  problem  should  be  conducted  in  order  to  establish 
the  magnitude  of  these  loads  and  establish  a  general  procedure  for  predict¬ 
ing  such  loads. 


1  odds  will  occur  in  regular  waves  in  the  tank  which  are  of  the  same  magnitude 
as  predicted  by  the  theory.  Also  we  need  to  determine  experimentally  the 
springing  loads  for  wave  groups  and  irregular  sea  waves.  The  springing  experi¬ 
ment  could  be  conducted  quite  inexpensively  by  using  a  vertical  cylinder  which 
can  pitch  about  a  fixed  horizontal  axis  at  the  pitch  center.  The  proper  re¬ 
storing  coefficient  pan  be  modelled  by  a  stiff  horizontal  spring  attached  at 
the  top.  Themost  important  aspect,  at  least  for  the  initial  experiment,  is  to 
model  correctly  the  physical  concept  of  the  springing  problem  and  to  do  it 
in  a  simple  way  so  that  the  correct  conditions  can  easily  be  modelled  and 
controlled. 

Furthermore,  improvements  are  needed  in  the  prediction  of  the 
extreme  values  for  both  the  surge  displacements  and  the  tendon  tension 
springing  loads.  Both  of  these  responses  are  natural  frequency  responses 
which  are  governed  by  nonl inear  viscous  damping  and  caused  by  nonlinear 
exciting  forces  with  frequencies  at  or  very  close  to  the  natural  frequency. 

Any  conventional  extreme  statistics  approach  as  used  for  linear  responses 
are  not  applicable. 

Therefore,  it  is  strongly  recommended  that  a  new  general  procedure 
should  be  developed  for  predicting  more  accurately  the  extreme  values  of  the 
slowly-varying  surge  motions  and  the  tendon  tension  springing  loads.  It  is 
anticipated  that  such  a  procedure  would  consist  of  the  two  following  parts: 

•  A  method  for  predicting  the  probability  and  the  statistics  for 
extreme  events  which  can  cause  linear  or  nonlinear  excitation  at 
or  near  prescribed  high  or  low  frequency  values. 

•  Nonlinear  methods  which  can  predict  with  sufficient  accuracy  the 
slowly-varying  horizontal  responses  and  the  high  frequency  spring¬ 
ing  loads  due  to  events  as  described  by  the  above-stated  method. 

It  is  believed  that  the  methods  presented  here  or  similar  methods  could  be 
used  as  a  starting  point  for  general  extreme  value  prediction  methods;  how¬ 
ever,  a  considerable  amount  of  improvements  and  extensions  would  be  required. 
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On  the  other  hand,  the  maximum  springing  loads  for  the  three  ten¬ 
dons  supporting  the  one  leg  may  for  natural  frequency  regular  wave  excita¬ 
tion  be  twice  as  large  as  it  is  for  a  TLP  with  all  sixteen  tendons  intact. 

1 . 4  Cone! ud ing  Remarks  and  Recommendations 

The  TLP  surge-motion  investigation  conducted  under  this  project 
seems  to  show  that  the  accuracy  of  surge  displacement  computations  could 
be  improved  by  a  better  viscous  force  model.  In  particular,  we  need  to 
improve  the  accuracy  of  the  prediction  of  the  nonlinear  viscous  forces  due 
to  interactions  between  the  different  modes  of  motion.  At  least,  improve¬ 
ments  are  needed  for  the  three  following  components: 

•  Additional  steady  surge  displacement  due  to  viscous  interactions 
between  the  steady  current  and  the  oscillatory  motions. 

•  Additional  viscous  damping  of  the  slowly-varying  natural  frequency 
surge  motions  due  to  viscous  interactions  between  these  motions  and 
the  steady  current  and  the  wave-frequency  motions. 

•  Steady  and  slowly-varying  nonlinear  viscous  "drift"  excitation  due 
to  viscous  interactions  between  the  different  wave  components. 

It  is  recommended  that  experiments  should  be  conducted  with  a 
single  large  vertical  cylinder  in  a  towing  tank.  Tests  should  be  conducted 
schematically  for  a  large  combination  of  current  speeds,  wave-frequency 
motions  and  slowly-varying  motions.  Most  importantly,  the  results  should 
be  made  generally  available. 

Similarly,  single  vertical  cylinder  experiments  are  also  recom¬ 
mended  for  the  springing  problem.  We  need  a  much  better  basic  understanding 
of  this  problem.  It  is  believed  that  a  single  cylinder  test  would  give  much 
more  useful  data  than  a  test  conducted  with  a  complete  TLP  model.  It  is 
extremely  important  to  determine  experimantal ly  if  tendon  tension  springing 


(c)  Sea-state  investigations 


Computations  performed  for  different  irregular  wave  condi¬ 
tions  seem  to  indicate  that  the  second-order  springing  loads  will 
be  substantially  larger  in  a  partially  developed  sea  state  than  in 
a  fully  developed  sea  state  with  the  same  significant  wave  height. 

(d)  Pitch  center  investigation 

The  pitch  center  investigation  results  show  that  a  very 
effective  way  to  reduce  the  second-order  tendon  tension  springing 
loads  is  to  move  the  pitch  center  closer  to  the  free  surface  level 
However,  the  study  shows  that  the  first-order  low-frequency  tendon 
loads  may  increase  by  moving  the  pitch  center  closer  to  the  free 
surface.  Therefore,  a  careful  investigation  of  the  first-  and 
second-order  tendon  load  results  is  required  before  any  judgment 
can  be  made  with  regard  to  the  pitch  center  location. 

(e)  TLP  size  investigation 

Springing  loads  were  predicted  both  for  the  2000-foot  TLP  #1 
design  and  for  the  3000-foot  TLP  #2  design.  A  comparison  of  these 
results  seems  to  indicate  that  the  springing  problem  is  much  more 
critical  for  a  deep-water  TLP  which  typically  has  a  higher  natural 
period. 

1.3.3  Results  for  Reduced  Number  of  Tendons 

The  investigations  presented  here  show  that  a  TLP  with  one  tendon 
removed  can  only  experience  an  external  load  50- 56 v  as  large  as  a  TLP  with 
all  of  the  tendons  intact  if  it  is  assumed  that  the  maximum  design  load  for 
a  tendon  is  twice  the  original  pretension  value. 

The  surge  motion  characteristics  for  a  TLP  with  one  less  tendon 
can  be  expected  to  be  quite  the  same  as  for  an  intact  TLP. 


Note  that  the  tendon  tension  load  due  to  first-order  regular  wave 
excitation  with  wave  frequency  equal  to  the  natural  frequency  is 
only  470  kips  per  tendon. 

If  TLP  #1  is  excited  by  two  simultaneous  regular  waves  with 
sum  frequency  precisely  equal  to  the  natural  frequency 

u),  +  =  to  =  2.508  rad/sec 

1  2  n 

the  additional  tendon  tension  load  due  to  second-order  wave  exci¬ 
tation  is 

(21 

t  -  780  kips  per  tendon 

Again  we  must  stress  that  these  values  must  only  be  considered 
as  upper  bound  estimates.  It  is  very  unlikely  that  such  precise 
regular  wave  events  will  ever  occur.  We  need  to  investigate  the 
probability  and  the  statistics  of  such  events  in  order  to  make 
better  estimates  of  the  extreme  springing  loads. 

(b)  Typical  operational  estimates 

For  a  typical  partially-developed  irregular  sea  condition 
with  significant  wave  height,  =  10',  the  significant  ampli¬ 
tude  of  the  wave-induced  nonlinear  tendon  tension  load  is 

t1/3  =  ^  kiPs  ^er  tend°n- 

This  value  has  been  obtained  by  a  new  irregular-wave  method  which 
includes  the  second-order  springing  excitation  due  to  the  inter¬ 
actions  between  the  different  wave  components.  It  is  believed 
that  this  method  will  give  reasonable  estimates  of  the  "average" 
tendon  tension  loads  for  typical  irregular-wave  conditions.  How¬ 
ever,  it  is  believed  that  extreme  values  may  not  be  accurately 
estimated  by  this  approach. 


Furthermore,  it  is  of  interest  to  note  that  the  predic¬ 
tions  show  that  if  wave  groups  with  perfect  tuning  do  occur 
in  the  ocean,  they  may  result  in  slowly-vary inq  surge  motions 
of  considerable  magnitude  and  probably  larger  than  those  due  to 
the  maximum  unsteady  wind  conditions  considered  here. 

1.3.2  Springing  Load  Results 

We  shall  here  refer  to  springing  loads  as  the  additional  tendon 
loads  caused  by  wave-induced  vertical  motions  at  or  near  the  natural  fre- 
quency  of  the  TLP.  The  high-frequency  vertical  motions  at  the  tendons  is 
due  to  natural-frequency  pitch  motions  about  the  virtual  mass  center.  The 
pitch  exciting  moment  is  caused  by  linear  and  nonlinear  wave-induced  hori¬ 
zontal  forces  acting  on  the  four  legs  at  a  location  close  to  the  undisturbed 
free-surface  level . 

Note  that  all  of  the  springing  load  results  presented  here  are 
for  zero  degree  heading.  It  follows  from  the  geometry  considerations  that 
it  can  be  expected  that  the  springing  loads  will  be  approximately  \fz  times 
larger  at  45  degree  heading  than  at  zero  degree  heading. 

(a)  Upper  bound  estimates 

For  TLP  #1  excited  by  regular  waves  with  period  precisely 
equal  to  twice  the  natural  period 

T  -  2T  =  5.012  sec, 
w  n 

the  additional  tendon  tension  load  due  to  natural  frequency  pitch 
motions  caused  by  second-order  wave-induced  excitation  is 

=  1020  kips  per  tendon. 


1-7 


fundamental  idea  is  to  use  conventional  finite  elements  to  approximate  the 
hydrodynamic  solution  in  a  fluid  region  near  the  body,  whereas  analytical 
representations  are  employed  outside  this  region.  By  using  a  variation 
formulation,  finite  elements  are  necessary  only  in  a  region  very  close  to 
the  body  resulting  in  significant  savings  in  computational  effort  (see 
Figure  2-1).  The  flexibility  and  versatility  of  the  finite  element  approach 
and  the  virtual  absence  of  tedious  analytical  preparation  make  the  HFEM  a 
viable  and,  in  many  applications,  superior  alternative  to  the  traditional 
integral  equation  method  using  Green  function  (source-distribution  method). 


Figure  2-1.  Typical  HFEM  grid  for  a  sample  TLP  configuration. 


The  various  theoretical  and  application  aspects  of  the  HFEM  are 
well  established  (Yue,  Chen  and  Mei,  1976;  Yue,  Chen  and  Mei ,  1978;  Aranha, 
Mei  and  Yue,  1979),  and  the  method  is  now  widely  accepted  among  offshore 
engineers.  Further  verification  of  the  improved  SAI  HFEM  code  is  also 
available  in  SAI  Capability  Report  -SAI -83-463-07 . 


2-9 


The  second-order  steady  surge  drift  force  is  formulated  using 
momentum  conservation  relationships  expressed  in  terms  of  the  far-field 
radiated  or  scattered  velocity  potential.  The  original  idea  is  due  to 
Maruo  (1960)  and  our  derivation  follows  closely  to  that  of  Newman  (1967). 

HFEM  Result 

The  HFEM  computer  program  has  been  used  for  all  of  the  results 
presented  in  this  section.  The  TLP  has  in  these  computations  been  repre¬ 
sented  by  the  four  vertical  legs  and  the  four  horizontal  pontoons.  Special 
attention  has  been  given  to  three-dimensional  effects  and  interactions 
between  the  major  members  (legs  and  pontoons).  Only  0°  heading  has  been 
considered.  All  of  the  HFEM  results  presented  in  this  section  are  for 
TLP  #1. 

Aside  from  the  four  vertical  legs  and  four  horizontal  pontoons, 
all  other  small  submerged  structural  members  and  risers  are  ignored  in  the 
diffraction  theory.  Such  members  constitute  only  about  5%  of  the  total  sub¬ 
merged  volume  and  their  typical  dimensions  are  so  small  compared  to  incident 
wave  lengths  that  diffraction  effects  can  be  ignored.  (A  term  corresponding 
to  the  hydrodynamic  added  mass  of  these  small  members  is,  however,  included 
in  the  equation  of  motion  as  a  small  correction.)  Since  the  vertical  dis¬ 
placement  of  the  TLP  is  much  smaller  than  that  in  surge,  its  contribution 
to  horizontal  drift  forces  is  unimportant  and  is  neglected  in  our  drift 
calculations. 

The  following  HFEM  results  for  TLP  #1  are  presented  as  function 
of  incident  wave  period: 

•  Figure  2-2a.  Added  mass  in  surge. 

•  Figure  2-2b.  Wave  radiation  damping  in  surge. 

•  Figure  2- 3a.  First-order  surge  exciting  force,  KAO. 


•  Figure  2-3b.  First-order  surge  exciting  force  amplitudes  for 

H/X  as  given  by  Figure  2-5. 

•  Figure  2-4a.  Second-order  drift  force,  RAO. 

•  Figure  2-4b.  Second-order  drift  force  amplitudes  for  H/X  as 

given  by  Figure  2-5. 


•  Figure  2-5.  Assumed  relationship  between  maximum  wave  steepness, 

H/X,  and  wave  period,  T. 

The  following  are  main  highlights  of  our  computer  results. 

•  For  moderately  long  waves  (T  ^  20  seconds),  the  surge  added  mass 
based  on  displaced  mass  (Morison's  formula  with  inertia  coeffi¬ 
cient,  Cm  =  2.0)  is  almost  30%  larger  than  predicted  by  the  HFEM 
code.  This  is  mainly  due  to  three-dimensional  and  interaction 
effects . 

•  For  exciting  force  and  linear  surge  motion,  most  of  the  effect 
of  multiple  legs  can  be  accounted  for  by  relative  phasing,  and 
the  effect  of  hydrodynamic  interaction  is  relatively  small.  The 
presence  of  the  pontoons  is  important  for  longer  waves,  T  >  8  sec. 
For  example,  for  T  =  14  sec,  the  pontoons  increase  the  exciting 
force  by  almost  20%. 

•  For  the  surge  drift  forces,  the  effect  of  hydrodynamic  interac¬ 
tion  is  significant,  and  the  results  based  on  a  single  leg  calcu¬ 
lation  are  inadequate.  The  drift  force  is  not  greatly  changed 

by  allowing  the  body  to  surge,  and  the  pontoons  have  almost  negli¬ 
gible  contribution.  The  heading  angle  of  the  TLP  is  found  to  have 
an  important  effect  on  the  drift  force,  but  this  aspect  has  not 
been  further  investigated  in  this  study. 
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figure  2-5.  Assumed  relationship  between  maximum  wave  steepness, 

H/A,  and  wave  period,  T. 

2 . 3  Ecju a tion  of  Motion 

A  computer  program  has  been  developed  which  predicts  the  nonlinear 
surge  motions  in  the  time  domain.  A  numerical  time- steppi ng  procedure  is 
used  in  which  the  time,  t,  is  advanced  in  small  time  steps,  At.  The  equa¬ 
tion  of  motion  which  is  solved  at  each  time  step  is 

(M  +  A)x  +  Bx  +  Cx  =  fy  +  +  f^.  (2-1) 

where  M  =  total  body  mass 
A  =  surge  added  mass 

B  -  surge  wave-radia tion  damping  coefficient 

C  =  nonlinear  surge  restoring  coefficient 

x  =  surge  acceleration 

x  =  surge  velocity 

x  =  surge  displacement 

fy  =  total  hydrodynamic  viscous  force 

fy  -  visions  wind-drag  force 

f^  -  total  wave- induced  exciting  force 
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The  total  mass  of  the  TLP  HI  design,  including  tendons  and  risers, 
is  (see  Appendix  A) 

M  =  2,600  *  103  slugs. 

The  linear  spring  constant,  the  value  of  C  as  x  0,  is 


=  18.32  kips/ft. 

The  added  mass.  A,  and  wave  damping  coefficient,  B,  have  been  com¬ 
puted  by  the  HFEM  code  (see  Section  2.2).  For  the  computations  presented 
here  the  infinite  period  values  (T  -*■  °°)  have  been  used: 

A(T  -*■«)  =  2,630  *  103  slugs 
and 

B(T  -*■«)  =  0. 


The  total  hydrodynamic  viscous  force,  f^ (t) ,  is  computed  at  each 
time  step  by  the  SAI  drag  formula 

V*)  =  Z^Pd{cds  lVsl  Vs^1+  f<VIVs*’  KC1)+  f(Vl  VSI  ’  KC2^ 

+  CD1  lVll  V10(y1>  KC1}  (2-2) 

+  CD2  |V2|  V2B(y2,  Kc2)  [1  +  f(Vj/V2,  Kcl)]j 


See  Salvesen,  et  al.  (1982)  for  more  details.  Note  that  the  drag  coeffi¬ 
cients,  Cq,  and  the  relative  velocity  components  vary  for  the  different 
sections  and  that  the  summation  is  over  all  of  the  sections  which  the  dif¬ 
ferent  members  are  divided  into.  For  comparison  purposes,  the  viscous  drag 
forces  are  also  computed  by  the  "quadratic  drag  formula", 


fV(t)  =  X^pD  c0  I V|  V 


2-16 


.  ..i-. . 
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where  V  is  the  fluid  relative  velocity  at  each  section.  The  above  formula 
is  often  referred  to  as  the  Mori  son  drag  force. 

The  viscous  wind-drag  force  is  predicted  by 

fw(t)  =  £%pCpA.  W|W|  (2-4) 

where  W  is  the  wind  speed.  The  maximum  steady  wind  speed  is 

Wq  =  70  knots, 

and  the  most  severe  unsteady  wind  is  assumed  to  be 

W(t)  =  70  +  7  •  cosui^t  (knots).  (2-5) 

Here  u>n  is  the  surge  natural  frequency  of  the  system.  The  maximum  steady 
wind  force  is 


fy  =  960  kips 


and  the  maximum  unsteady  wind  force  is 


fw(t)  =  960  kips  +  190  kips  •  cosu)nt. 


(2-6) 


The  total  wave-induced  exciting  force  is 


fE(t)  =  Fo  +  EFi  cos  (to  - 1  +  c.j )  +  £Fj  cos(Ujt+ej)  (2-7) 


where 


=  steady  wave  ^jrift  force 


=  amplitude  of  the  primary  wave  exciting  force  components  which 
has  the  same  frequency  as  the  frequency  of  wave  components 


2 . 4  Regular-Wave  Surge  Displacement  Results  for  TLP  fll 

The  surge  displacement  in  regular  waves  as  computed  here  will 
finally  come  to  a  steady-state  condition  after  the  initial  transient  motions 
have  died  out.  The  final  steady-state  surge  motions  x(t)  in  regular  waves 
may  be  expressed  as 

x(t)  =  XQ  +  Xj  cos^t  +  X2  cos^t  (2-10) 

where  XQ  =  mean  displacement  due  to  wind,  wave  drift  and  current 
(including  body  and  wave  velocity  interactions) 
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Xj  =  amplitude  of  the  frequency  components  of  the  primary  first- 
order  motions  due  to  primary  wave  excitation  only  (primary 
motions  are  dominated  by  inertia  effects  and  are  not  affected 
by  current,  wind  or  the  slowly-varying  motions) 

=  amplitude  of  the  slowly-varying  motion  components  due  to 
unsteady  wind  (its  magnitude  is  affected  by  current  and  the 
primary  motions  through  the  viscous  damping  force). 

The  mean  displacement  Xo  and  the  amplitude  of  the  slowly-varyinq 
motions  X^  shall  be  further  divided  into  parts  which  describe  tne  separate 
conti rbutions  due  to  wind,  wave  drift,  current  and  current-motion  interac¬ 
tions.  The  mean  surge  displacement  can  be  divided  into  the  following  five 
components: 


v  VW  A  VD  x  VC  x  VCP  .  VCS 
Xo  =  Xo  +  Xo  +  Xo  +  Xo  +  Xo 


(2-11) 


where 


XQ  =  displacement  due  to  steady  wind 


X^  =  displacement  due  to  steady  wave-drift 


X^  =  displacement  due  to  current  alone 
CP 

XQ  =  displacement  due  to  steady  viscous  force  caused  by  inter¬ 
action  oetween  current,  primary  body  and  wave  particle 
velocities 


CS 

XQ  =  additional  displacement  due  to  steady  viscous  force  caused 
by  interactions  between  current  and  the  slowly-varying  body 
motions. 


The  amplitude  of  the  slowly-varying  motions,  X^,  can  be  divided  into  the  fol 
lowing  three  components 
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where  =  amplitude  of  the  slowly-varying  motions  due  to  unsteady 

wind  with  no  current  present  and  no  primary  motions 

IC 

X2  =  increase  (or  decrease)  in  amplitude  due  to  current 
I P 

=  additional  increase  (or  decrease)  in  amplitude  due  to  the 
primary  motions. 

An  investigation  of  these  different  surge  displacement  components 
is  presented  in  Appendix  B.  The  surge  displacements  are  computed  using  both 
the  SAI  and  the  quadratic  drag  formula  for  the  0°  heading  case.  The  TLP  is 
excited  by  maximum  current,  maximum  regular  waves  (H  =  80.0  feet  and  T  =  14.0 
seconds)  and  unsteady  wind  (W  =  70.0  +  7.0  cosajnt  (knots)  where  wn  is  the 
natural  surge  frequency). 

The  computed  time  domain  surge  displacement  for  this  wind,  wave 
and  current  condition  is  presented  in  Figure  2-6.  Results  are  shown  using 
both  the  SAI  drag  formula  and  the  quadratic  formula.  It  is  seen  that  the 
SAI  drag  formula  predicts  that  the  three  main  surge  displacement  components 
in  this  case  are 


Steady  surge 

Amplitude  of  primary  motion 
Amplitude  of  slowly-varying  motion 


X  =  70.5' 

o 

Xj  =  22.1' 

X2  =  9.7' 


Total  maximum  surge  displacement 


max 


102.4' 


When  the  quadratic  formula  is  used  the  components  are 
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Figure  2-6.  Computed  nonlinear  time-domain  responses  for 
TLP  excited  by  current,  unsteady  wind  and 
maximum  design  wave  with  0°  heading. 


Steady  surge  XQ  =  96.7' 

Amplitude  of  primary  motion  Xj  =  22.1' 

Amplitude  of  slowly-varying  motion  X2  =  8.9' 

Total  maximum  surge  displacement  Xmax  =  127.7' 

We  shall  now  present  results  for  the  additional  individual  com¬ 
ponents.  The  different  surge  displacement  components  have  been  determined 
by  performing  runs  with  the  time  domain  surge  motion  computer  program  using 
different  combinations  of  the  above-stated  current,  wave  and  wind  conditions. 
The  results  for  each  of  these  computer  runs  are  presented  in  Appendix  B.  A 
summary  of  the  results  is  given  in  Table  2-2. 

The  most  important  final  conclusions  which  can  be  drawn  from  the 
results  presented  in  Table  2-2  and  in  Appendix  B  are: 

t  Steady  surge  displacement  is  14.4  feet  when  the  TLP  is  excited 
by  current  alone.  When  the  current  and  the  maximum  design  wave 
are  both  present  the  quadratic  drag  formula  predicts  an  increase 
of  the  steady  displacement  by  27.9  feet  due  to  interactions  be¬ 
tween  the  current,  the  wave  particle  motions  and  the  body  motions. 
Results  obtained  by  the  SAI  formula  show  only  a  1.8-foot  increase. 
Furthermore,  the  results  in  Table  2-2  show  that  if  slowly-varying 
motions  are  also  present,  then  the  steady  surge  is  unchanged  when 
the  quadratic  formula  is  used  and  decreased  by  0.1  feet  when  the 
SAI  formula  is  used.  Note  that  there  is  an  additional  steady 
surge  component  resulting  from  integration  of  the  viscous  forces 
over  the  instantaneous  wetted  surface.  This  component  which  is 
equal  to  about  3-4  feet  is  not  included  in  Table  2-2. 

•  Amplitude  of  the  primary  wave-induced  surge  motions  is  22.1  feet. 
Computations  have  shown  that  the  wave-frequency  motion  is  affected 
only  slightly  by  viscous  drag  forces.  These  motions  are  essen¬ 
tially  completely  dominated  by  inertia  effects.  This  implies  that 
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f i cj u r e  2-8.  Damping  factor,  c,,  as  function  of  surge  exciting 

force  amplitude  and  as  function  of  surge  displace¬ 
ment  amplitude,  X^. 

2./  TIP  Size  Investigation 

In  this  section  we  shall  compare  surge-motion  results  for  the  TLP 
■‘I  on c  7l?  -2  designs.  Detailed  descriptions  of  the  two  TIPs  are  included 
:r.  creix  A.  Dote  that  TLP  n?_  is  1.2  times  larger  than  TLP  fl,  and  that 
ILL’  -I  is  assumed  to  be  designed  for  2000-foot  depth  whereas  TLP  ?2  is 
•.esi^ned  for  3000-foot  depth. 


Tne  same  iiiTM  potential-flow  data  as  presented  in  Section  2.2  and 
tie  ti:;  c- domain  regular  wave  results  as  presented  in  Section  2.4  and  Appendix 
R  obtained  for  TLP  *1  have  also  been  obtained  for  TLP  f 2.  i-iowever,  here 
we  shall  not  present  all  of  the  detailed  results  for  TLP  v2.  Wc  shall  only 


•  sec. 


/If)  kij 


Now  using  f  igure  2-1'.,  «-fjua  I.  ion  (2-35)  is  easily  solved  by  iteration  <1  ru  1  the 
results  are 

-  2.6  feet  and  l  ~  0.06. 

1/  J 

The  final  surge  displacement  results  for  the  sample  irregular 
sea  condition  are 

Steady  Displacement 

Due  to  wind 

Due  to  wave  drift 

Primary  wave-frequency  motions 

significant  amplitude, 

S 1  owl y- vary ing  natural-frequency  motions 

(?) 

Significant  amplitude, 

Total  "significant"  displacement  11.6' 

If  we  r.ac  included  mu x  imum  current,  this  would  have  resulted  in  an  additional 
•  ‘t.vdy  surge  displacement  of  about  12  feet  ano  probably  have  reduced  the 
y- varying  motion  amplitude  by  a  factor  of  0.5.  Hence,  in  the  same 
irregu \uc  sea  condition  with  current,  the  total  "significant"  surge  displace- 
'•nt  will  do  approximately  2d  feet. 


5 .  >' ; 1 
2.3' 

0 . 9 1 

2.0‘ 
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where  A  =  0.0081  y  and  R  ~  33.56  h^-  The  significant  wave  height, 

10  feet.  The  average  mean  period  for  this  is  Tj  -  6.73  seconds.  We  shall 
assume  that  a  typical  wind  speed  for  this  sea  condition  is 


W  =  39  ft/sec. 

For  this  sample  case  we  shall  assume  that  the  current  is  zero. 

By  equations  (2-19),  (2-20)  and  (2-21)  and  numerical  HFEM  results 
given  in  Figures  2-2  through  2-5,  we  have  that  the  significant  first-order 
surge  motion  amplitude  is 


X(1)  -  2  0  o(1) 
l/3  °x 


0.92  feet. 


(2-33) 


The  steady  drift  force  is  by  equation  (2-22) 


F(2)  =  41  kips, 
o 


By  performing  computations  with  the  nonlinear  surge  time-domain  code  repre¬ 
senting  the  wave-frequency  surge  motion  by 


x^(t)  =  0.92’  •  cos t 


(2-34) 


we  nave  ootained  the  nonlinear  relationship  between  the  damping  and  the 
amplitude  of  the  natural-frequency  slowly-varyi ng  surge  motions.  This 
relationship  is  shown  graphically  in  Figure  2-8. 


From  equation  (2-28)  it  follows  that  the  significant  amplitude 


of  the  slowly-varyi  ng  motion,  for  an  irregular  sea  is 


I2  •  x$)  =  42><“'*>  Ht-3 


(2-35) 


where  by  applying  equation  (2-25)  we  have  that 


w”  •  •  r  <.r,('  fir, l  f  r.omponon  t.  is  o  I  imvir  wave  •  f  req,.ency  component  and 

‘.hr‘  second  component.  is  a  s lowly-varying  second-order  component  with  fre¬ 
quency  equal  to  the  natural  surge  frequency. 

We  snail  assume  that  for  an  irregular  seaway,  an  equivalent 
linear  viscous  damping  factor  can  be  used  and  that  this  damping  factor  is 
tne  same  us  the  damping  factor  for  natural  frequency  surge  motions  with 
amplitude  equal  to  the  significant  second-order  surge  motion,  i.e. 

x(2)(t)  =  coswnt. 

Tr;e  nonlinear  relationship  between  the  significant  surge  di  splaccr.ient  anipli- 
f  ?  i 

tnde,  X,,,,  and  the  damping  factor,  q,  can  then  be  obtained  by  running  the 
1/  -J 

time-domain  surge  motion  computer  program  for  several  exciting  force  values. 

However,  since  the  viscous  damping  depends  on  nonlinear  interac- 
Ti.v,  resulting  from  current  and  wave-frequency  motion,  the  proper  current 
and  wave-frequency  motions  must  be  included  in  the  time-domain  computations. 
Wo  snail,  for  simplicity,  assume  that  the  wave-frequency  motions  effect  on 
tne  viscous  damping  can  be  represented  by  using  a  single  wave-frequency  com¬ 
ponent  given  by 

xU'(t)  --  xjj]  coswjt  (2-31) 

f  1  ' 

.-.ser,?  ....  is  the  average  mean  frequency  for  the  given  sea  spectrum  and 
,t  tne  significant  surge  amplitude  for  the  first-order  motion. 

2.6.2  Sample  Results 

We  shall  assume  that  a  typical  seaway  is  represented  by  a  fully 
d.'V'  I--:.,;.'.  Pierson-Moskowi  1/  spectrum 


c,  =  B/2Au)n.  (2-29) 

Note  that  for  surge  motions  at  and  near  the  natural  frequency,  the  damping 
is  nonlinear  and  depends  on  the  magnitude  of  the  surge  motions.  Furthermore, 
the  surge  damping  depends  on  complicated  nonlinear  viscous  interaction 
between  the  si owly-varying  motions  and  the  steady  current  and  also  the  £.i_TJmt- 
order  wave-frequency  motions . 

We  shall  here  obtain  the  nonlinear  relationship  between  the  damp¬ 
ing  factor,  and  the  surge  motion  amplitude  by  using  our  nonlinear  time- 
domain  surge  motion  computer  program.  This  computer  program  assumes  that 
the  surge  motions  only  have  two  frequency  components, 

x(t)  =  Xq  +  X^cosujt  +  X^cosunt  (2-30) 
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Note  that  the  o  elded  mu  v»  and  damp  ;  r:(|  coef  f  ic  i  ent s ,  A  arid  b ,  are  f  unctions  of 
frequency,  <a,  as  shown  in  Figures  2-2a  and  2-2b.  I  no  variance  (:.e..  the 
square  of  the  RMS  value,  a)  of  the  first-order  sureje  displacement  is  given 

by 


d/^(u))  due 
X 


(2-21) 


The  mean  steady  second  order  drift  force  is 


p  ( 2 )  _  9 

o 


rx 

J  F^(w)  •  4a(w)  dw 


(2-22) 


where  Fu^U)  is  the  steady  second-order  drift  force,  RAO,  for  regular 
(see  figure  2-4a). 


waves 


"lie  surge  force  power  spectrum  due  to  second -order  wave- induced 


excitation  is 


u> 

'/*a{w  + ! 


)  *>  - 1 


r)  |»U  +  ,  w  -  u‘-  )\  duj.  (2-23) 


here  ,  v^)  is  the  second-order  drift  force  RAO  for  excitation  by  two 

simultaneous  waves  with  frequencies,  <Uj,  and  If  <a*  is  small,  we  may  apply 

Newman's  approximation  and  assume  that 


£(&:  *  -  (2‘*)  =  F(Z)(U>) 


,  2-24) 


d  n  r*. ,  ^  ..j  r  1 1  *  j  Pi  i  or  t.j  ,  t,  r  d  t 


’}(m)  =  2  /(»)  F(2)(oi)]  do). 


(2-25) 


•,i.n;o  power  sp.e«- tren  due  to  second-order  excitation  is 


'?),  s  A2),  v  :  l:/  V  2 


Ol!/ 


(2-26) 
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X.,  40.0  fort  . 


Hence  if  wave  groups  with  perfect  tuning  do  occur  in  the  ocean,  they  wni 
result  in  slowly-varyi ng  surge  motions  which  are  considerably  larger  than 
those  due  to  the  maximum  unsteady  wind  condition  considered  here. 

2.6  I r r egular  Wave  Results  for  TLP  // 1 

The  main  objectives  with  the  irregular  wave  investigation  are  to 
present  a  method  for  predicting  the  "average"  surge  motions  in  an  irregular 
seaway  and  to  present  results  for  a  typical  sea  condition.  This  method  is 
not  recommended  for  predicting  extreme  surge  responses.  It  is  believed  that 
the  extreme  responses  can  be  predicted  more  accurately  by  computing  the 
nonlinear  time-domain  responses  for  particular  extreme  events  which  are 
expected  to  result  in  large  surge  responses. 

2.6.1  Theoretical  Formulation 

We  shall  assume  that  the  seaway  is  unidirectional  and  represented 
by  a  known  power  spectrum 

(w ) 

a v 

The  surge  displacement  power  spectrum  due  to  first-order  wave-induced  exci¬ 
tation  is 


41}(J  =  <DaU)  jF(1)(u))|2  •  J H ( o>)  | 2  (2-19) 

where  F^U)  is  the  first-order  surge  exciting  force,  RAO  (see 
Figure  2-3a) 

H(m)  is  the  "linear"  surge  displacement,  RAO  given  by 

1 

HU)  = - ? -  (2-20) 

(C  -  Aw  )  +  iBw 
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U) 


•,(•<. orni  unirr  ■ . !  < j  f  1  y  drift  fnri  <  lor  risi.i.ir  waves  ol 
ur.it.  wave  amplitude 


F^'  -  amp!  i  t  ude  of  si  owly- vary  i  ncj  force  for  two  wave  components 


M2) 


amplitude  of  t.he  "springing"  force  for  two  wave  components 


The  results  presented  in  Section  2.2  show  that  the  maximum  drift 


force  RAO  is 


(2) 


2.5  kips  when  T  -  9.3  sec. 


donee  we  have  that  the  second-order  excitation  for  the  wave-grouping  condi¬ 
tion  considered  here  is 


f^(t)  =  390  kips  i  390  kips  •  cosu.  t. 


inis  exciting  force  will  result  in  a  steady  surge  displacement 


390 

X..  -  =21.3  feet. 


o  18.32 


(2-18) 


and  a  slow iy- varying  surge  with  an  amplitude  (estimated  from  the  results 
presented-  in  Appendix  B) 


~  58  feet. 

‘ n i is  me  estimated  amplitude  when  there  is  no  current  and  when  damping 
u.i.  -u  interactions  witn  the  primary  motions  is  ignored.  If  current  '  s 
present  tr<  ar..pl i tucie  will  be  decreased  due  to  additional  viscous  damping 
forces  caused  by  interactions  between  trie  current  and  the  oscillatory 

;  (.  1  Of  i f-  . 


it  tin  ,/c  rec.lisc  l  nut  the  maximum  amplitude  for  the  slowly- 
otio-r.  resul  ting  ."ro.;-.  tnc  maximum  unsteady  wind  eur.u  i  Lion  is 
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Correct  to  the  second  order  in  wave  amplitude,  the  exciting  force 
due  to  two  regular  waves  is 


f(t)  =  f(1)(t)  <  f(2)(t)  +  .  .  . 


(2-13) 


where  the  first-order  force 


f^(t)  =  F^Uj  coswjt  +  F^^a2  cosuj^t 


(2-14) 


and  the  second-order  force  is 


f(2)(t) 


=  F{2)(Wl)  •  a2  +  F(2)(co2)  • 


_  {2) 

+  2a^a2  F'  w2)  •  cos (w2  -  co^)t 

4  (  2  ) 

+  2a jU2  Fv  ^  (w^ »  (i>2 )  •  cos(io1  +  <»>2)t 


(2-15) 


+  a2  V2*^)  cos2ujjt  +  a £  +  F^2^(u>2)  cos2w2t 


We  are  here  only  interested  in  the  steady  and  the  slowly-varying  part  of  tne 
second-order  force  which  may  be  written  as 


A  2)m  _  A  2)/  >  2  ,  p(2),  >  2 

t  (t)  -  t-  v^j)  •  +  F  *  «2 


(?)  W1  +  w2 
+  2a^a2  Fv  '  (■  ■  -  2~  )  •  cos(to2  -  w^)t 


(2-16) 


wr.ere  the  Newman  approximation  has  been  applied.  The  Newman  approximation 
states  that 


■f'%,  „2)  -  F (2)("1^>  ♦  0(„2  - 


(2-17) 


Here 


amplitude  of  first-order  exciting  force  per  unit  wave 
amp! i tude 


2.(j  bu r< n:  Mol. ion'.  Due  to  Wave  fit'ou | » i h«j 

We  shall  consider  excitation  due  to  "wave  groupi tn;"  wim  perfect 
tuning.  That  is  two  sin.ul taneous  regular  waves  with  difference  frequency 
equal  to  the  surge  natural  frequency  and  with  ii,oan  frequency  equal  to  the 
frequency  for  which  tne  steady  surge  drift  force  is  maximum.  This  particu¬ 
lar  wave  condition  will  result  in  the  largest  possible  s lowly-varying  wave- 
induced  excitation  since  all  of  the  wave  energy  is  concentrated  at  the  maxi¬ 
mum  wave-drift  frequency.  However,  it  must  be  recognized  that  it  has  not 
been  established  with  any  degree  of  confidence  that  wave  groups  with  perfect 
tuning  will  actually  occur  in  the  ocean.  The  wave-grouping  results  pre¬ 
sented  here  should  be  considered,  therefore,  as  an  estimate  of  the  upper 
limit  for  the  wave-induced  slowly-varying  excitation  force. 

We  shall  consider  two  waves: 


ies  that 


o  .  aS  f  i  e  t. . 


Tattle  2-  3 .  '.urge  d  i |  >  I  «>  c  cun  ri  t Ini  two  i  ■ »  j  i  «  -jih  •  inquiur  wave 
condition',  for  0°  heading. 


Maximum  Design 
Wave 

H  =  30.0  ft 

T  =  14.0  sec 

Maximum  Drift 
Wave 

H  =  46.4  ft 

T  =  9.5  sec 

SA1 

.  Draa _ 

Quad 

Drag 

SAI 

Oraq 

Quad 

Drag 

r 

Steady  Current  Displacement,  XQ 

16.1 

42.3 

15.2 

35.5 

w 

Steady  Wind  Displacement,  XQ 

52.5 

52.5 

52.5 

52.5 

Steady  Drift  Displacement,  X^ 

2.0 

,0 

64.1 

64.1 

Wave  frequency  Amplitude, 

22.1 

22.1 

1.7 

1.7 

Slowly- Varying  Motion  Amplitude,  X^ 

9.7 

8.9 

18.8 

12.8 

Total  Maximum  Surge  Displacement 

102.4 

127.7 

152.3 

166.6 

The  results  summarized  in  Table  2-3  show  that  the  largest  surge 
displacement  will  occur  for  the  0°  heading  angle  with  "maximum  drift  wave" 
and  maximum  current.  When  the  quadratic  drag  formula  is  used,  the  maximum 
displacement  is 


X  a  -  166.6  feet, 
ma  x 

Note  that  the  results  presented  in  Table  2-3  show  that  for  the  maximum 
design  wave  case  the  surge  displacement  (using  quadratic  drag)  is 

Xmax  =  127 • 7  feet 

which  is  substantial ly  smaller  than  the  surge  displacement  for  the  "maximum 
drift  wave"  case. 
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Figure  2-7.  Surge  displacements  for  "maximum  drift  wave"  condition 
with  0°  heading  (Condition  2). 
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The  results  presented  above  indicate  that  the  SAI  drag  formula 
predicts  values  of  viscous  forces  due  to  interactions  between  current  and 
wave-frequency  motions  that  are  much  smaller  than  the  values  predicted  by 
the  quadratic  viscous  force  formula.  They  also  show  that  the  total  surge 
motions  predicted  by  the  SAI  drag  formula  will,  in  general,  have  smaller 
steady  displacements  than  the  motions  predicted  by  the  quadratic  drag  formula. 

It  is  difficult  to  determine  the  particular  environmental  condi¬ 
tion  and  heading  angle  which  will  result  in  the  largest  maximum  surge  dis¬ 
placement.  Note,  for  example,  that  the  current  may  decrease  the  amplitude 
of  the  slowly- varying  motions  by  a  magnitude  larger  than  the  steady  dis¬ 
placement  due  to  current.  In  such  a  case,  the  current  will  actually  result 
in  a  decrease  in  the  total  surge  displacement.  Several  computer  runs  have 
been  performed  in  order  to  establish  the  most  extreme  condition.  Of  par¬ 
ticular  interest  are  the  following  two  regular  wave  conditions  for  which 
results  are  summarized  in  Table  2-3. 

Condition  1 

0°  heading,  maximum  regular  wave  with  H  =  80.0  feet  and  T  =  14.0  sec¬ 
onds,  maximum  unsteady  wind  and  maximum  current.  Motions  computed  by 
SAI  and  quadratic  drag.  Results  for  the  surge  motions  as  functions  of 
time  are  presented  in  Figure  2-6. 

Condition  2 

0°  heading,  "maximum  drift  wave"  with  H  =  46.4  feet  and  T  =  9.5  sec¬ 
onds  and  maximum  unsteady  wind.  Surge  motion  plots  for  this  condition 
are  presented  in  Figure  2-7. 

The  "maximum  drift  wave"  in  Condition  2  refers  to  the  regular 
wave  which  will  produce  the  largest  steady  wave-induced  drift  force  for  the 
heading  in  question.  It  is  assumed  in  these  cases  that  the  wave  amplitude 
is  determined  from  the  maximum  wave  steepness  given  by  H/A  =  0.10.  Note 
that  the  maximum  theoretical  wave  steepness  for  nonbreaking  waves  is  H/X  = 

0.14  and  that  waves  with  steepness  H/X  >  0.12  have  in  many  cases  been 
recorded . 

*  As  shown  in  Figure  2-4b,  the  maximum  draft  force  in  regular  wave  is  1,200 
kips  and  occurs  at  T  =  9.5  sec. 


the  wave-frequency  motions  can  be  predicted  very  accurately  by 
linear  theory  ignoring  viscous  damping  and  interactions  due  to 
current  and  slowly-varying  motions. 


•  Amplitude  of  the  slowly-varying  motions  due  to  unsteady  wind  exci¬ 
tation  alone  40.9  feet.  If  current  is  present  in  addition  to 
the  unsteady  wind,  then  the  amplitude  of  the  slowly- varyi ng 
motions  is  predicted  to  decrease  by  lb. 4  feet  when  the  SAI  drag 
formula  is  used  and  25.4  feet  when  the  quadratic  formula  is  used. 
If  wave-frequency  motions  are  also  present,  then  the  SAI  formula 
predicts  an  additional  decrease  of  15.8  feet,  whereas  the  quadrati 
formula  predicts  an  additional  6.6-foot  decrease.  The  total 
decrease  in  the  slowly-varying  amplitude  due  to  current  and  wave 
frequency  motion  is  31.2  feet  when  the  SAI  drag  formula  is  used 
and  32.0  feet  when  the  quadratic  formula  is  used. 

Table  2-2.  Surge  displacement  components  for  TLP#1  excited  by  current, 
maximum  design  regular  waves  and  unsteady  wind  with  0° 
heading  angle  (values  given  in  feet). 


SAI 

Drag 

formula 

Morljon 

Drag 

formula 

SAI 

Drag 

formula 

MoHson 

Drag 

formula 

5A1 

Drag 

Formula 

Mori  son 
Drag 
formula 

Current  alone 

ia.a 

mm 

Due  to  current 
t  Interactions 

16.1 

42.3 

Increase  due  to 
wave- frequency 
motion 

1.8 

27.9 

Steady 

1  displacement 

70.5 

96.7 

Increase  due  to 

Slowly-varying 

motion 

■ 

Due  to  wind 

62.5 

52  5 

_ 

2.0 

MM 

Wave-f rpouency 
ampl Itude 

22.1 

22.1 

Slowly-varying  * 
excitation  alone 

40.9 

40.9 

Slowly-varying 
ampl 1 tude 

9.7 

8.9 

Increase  due  to 
current 

-15.4 

-25.4 

Increase  due  to 

wave-frequency 

motion 

-15.8 

-6.6 

Total  surge 
displacement 

102.4 

127.7 

Slowly-varying  due  to  unsteady  wind  excitation  alone. 


summarize  the  Til1  !!?  surge  motion  results  and  nuke  comparisons  helween  the 
two  platforms.  Note  that  no  irregular  wave  resuli-,  have  hern  obtained  for 
TLP  112. 

Table  2-4  shows  the  main  surge  motion  component  for  TLP  H2  excited 
by  the  maximum  design  wave  and  the  regular  wave  condition  which  results  in 
the  maximum  steady  drift  exciting  force  (referred  to  as  the  "maximum  drift 
wave").  Note  that  the  maximum  design  wave,  the  maximum  current  and  the 
maximum  wind  conditions  are  assumed  the  same  for  TLPs  ft  1  and  f>2.  Comparing 
tne  results  presented  in  Table  2-4  with  the  TLP  #1  results  shown  in  Table 
2-3,  it  is  seen  that  the  surge  motions  for  TLP  tt2  are  slightly  larger  than 
for  TLP  ft  1.  However,  if  we  normalize  the  surge  motions  with  respect  to 
water  depth,  the  TLP  ft 2  motions  are  in  general  smaller  than  for  TLP  ft  1.  It 
is  not  possible  to  make  any  specific  comparisons  between  the  the  two  TLPs. 

We  have  here  assumed  that  the  dimensions  for  TLP  il 2  are  1.2  times  the 
dimensions  for  TLP  #1.  However,  it  is  difficult  to  evaluate  if  1.2  is  the 
correct  scaling  between  a  platform  operating  in  2000-foot  depth  and  a  plat¬ 
form  operating  in  3000-foot  depth.  Therefore,  we  must  be  careful  in  drawing 
any  general  conclusions  from  a  comparison  between  the  surge  motions  for 
these  two  platforms. 

Table  2-5  shows  the  surge  displacement  components  for  TLP  “2 
excited  oy  current,  maximum  design  wave  and  unsteady  wind.  Comparing  these 
results  with  the  results  presented  in  Table  2-2  for  TLP  ttl  it  is  seen  that 
tne  viscous  interaction  effects  as  computed  by  the  Morison  quadratic  drag 
ror:,'u  i  a  and  the  SAI  drag  formula  are  about  of  the  same  relative  magnitude, 
.uwever,  wo  do  note  that  the  additional  damping  of  the  slowly-varying  motions 
dw  to  viscous  interaction  effects  between  the  slowly-varying  motions  and 
the  current  are  substantially  larger  for  TLP  It 2  than  for  TLP  1 .  When  the 
!>A»  drag  formula  is  used  this  additional  damping  effect  is  2.5  times  larger 
for  TLP  ■■■2  than  for  TLP  ?i,  and  when  the  Morison  drag  formula  is  used,  it 
is  i.o  ir.es  larger,  it  is  difficult  to  pinpoint  any  specific  reason  for 
tnese  differences. 
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SPRINGING  TENDON  LOAD  INVESTIGATION* 


3. 1  The  TLP  Springing  Problem 

3.1.1  General  Discussion  of  the  Springing  Problem 

Linear  and  nonlinear  wave-induced  high-frequency  springing  excit¬ 
ing  forces  can  result  in  large  oscillatory  loads  in  the  tendons  at  the  nat¬ 
ural  frequency  of  the  TLP.  The  natural  periods  for  heave  and  pitch  motions 
for  most  deep-water  TLPs  are  typically  between  2  to  3  seconds.  For  the 
discussion  in  this  section  let  us  assume  that  the  vertical  motions  natural 
period,  T  ,  equals  2.5  seconds. 

Regular  wave  excitation  shal  1  be  considered  first.  The  wave- 
induced  springing  exciting  forces  can  be  divided  into 

•  first-order  forces  resulting  from  a  wave  with  period,  Tw,  equal 
to  the  natural  period  of  the  TLP: 

T„  =  T  =  2.5  seconds 
w  n 

•  second-order  forces  resulting  from  a  wave  with  period  equal  to 
twice  the  natural  period  of  the  TLP: 

T,  =  2T  =5.0  seconds 
w  n 

•  third-order  forces  resulting  from  a  wave  with  period  equal  to 
three  times  the  natural  period  of  the  TLP: 

T  =  3T  =  7.5  seconds 
w  n 

anc  so  an  for  higher  order  excitations. 

*  All  of  the  springing  loan  results  presented  in  this  section  are  for  *ero 
,-;r.-e  heading.  The  springing  loads  for  45°  heading  can  be  expected  to 
,  •  J'i.  times  larger. 
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It  is  important  to  recognize  that  the  second-  and  third-order 
excitations  are  caused  by  waves  which  have  a  larger  period  and  consequently 
are  much  longer  and  have  much  more  energy  than  the  waves  giving  rise  to 
the  linear  excitation.  If  we  consider  a  constant  wave-height/wave-length 
ration,  H/X,  it  follows  that  the  waves  causing  the  second-  and  third-order 
excitation  have  _16  and  81  times,  respectively,  the  energy  of  the  waves 
causing  the  linear  excitation.  It  can  be  expected,  therefore,  that  the 
second-  and  third-order  springing  forces  may  be  much  larger  than  the  first- 
order  forces. 

Another  important  fact  is  that  the  springing  loads  in  the  tendons 
will  mainly  be  caused  by  pitching  motions  resulting  from  horizontal  wave 
excitation.  The  vertical  springing  excitation  due  to  short  waves  is  so 
small  that  it  can  practically  be  ignored. 

For  a  typical  2000-foot  depth  TLP  configuration  excited  by  regular 
periodic  waves  with  H/X  =  0.075,  our  computations  show  that  it  can  be  ex¬ 
pected  that  the  first-order  surge  springing  exciting  force  is  of  the  order 
of  magnitude  of 


fjlj  =  50  kips 

and  that  the  second-order  springing  force  may  be  of  the  order  of  magnitude  of 

=  200  kips 

and  that  both  act  very  close  to  the  free  surface  (about  three  feet  below). 

The  effective  pitch  center  for  a  TLP  can  be  expected  to  be  about  20  feet 
below  the  free  surface  resulting  in  a  pitch  moment  arm  of  about  17  feet. 
Hence,  the  amplitudes  of  the  first-  and  second-order  pitch  springing  moments 
may  be  of  the  order  of  magnitude  of 
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fr  =  800  kips -feet 
o 

=  3400  kips-feet 

Note  that  the  waves  causing  the  first-order  excitation  have  a  period  equal 
to  the  natural  period,  whereas  the  waves  causing  the  second-order  excita¬ 
tion  have  a  period  equal  to  twice  the  natural  period  of  the  TLP. 

Another  critical  aspect  is  that  there  is  very  little  viscous  damp¬ 
ing  and  practically  no  wave  damping  for  the  high-frequency  pitch  motions. 

The  amplification  factor  is  very  large  resulting  in  tendon  loads  which  may 
be  of  the  order  of  magnitude  of  600  kips  per  tendon.  Figure  1-3  shows 
scnematical ly  the  total  springing  surge  exciting  force  of  250  kips  and  the 
resulting  5,000  kips  load  on  the  eight  left-hand  side  tendons  and  -5,000 
kips  load  on  the  eight  right-hand  side  tendons.  These  are  approximate  results 
for  a  2000-foot  depth  TLP  with  typical  dimensions  as  shown  in  the  figure. 


1 

:j,!»  ;12:*250  kips 

Total  Springing 
inciting  Force 


We  have  at  the  present  time  made  no  attempt  to  predict  the  third- 
or  higher-order  springing  excitation.  It  is  believed  that  the  higher-order 
excitations  may  be  as  important  as  the  second-order;  however,  we  need  to 
have  a  better  understanding  of  the  second-order  excitation  and,  in  particu¬ 
lar,  as  it  applies  to  realistic  irregular  sea  conditions  before  we  can 
advance  to  the  more  complicated  higher-order  excitation  problem. 

3.1.2  Problem  Formulation 

We  shall  here  refer  to  springing  loads  as  the  additional  tendon 
loads  caused  by  wave-induced  vertical  motions  at  or  near  the  natural  fre¬ 
quency  of  the  TLP.  Since  the  wave  lengths  for  both  the  first-  and  second- 
order  springing  wave  excitations  are  smaller  or  of  the  same  order  of  magni¬ 
tude  as  the  depth  of  the  TLP,  it  can  be  concluded  that  the  heave  springing 
excitation  forces  are  small  and  can  be  ignored.  The  springing  loads  in  the 
tendons  at  or  near  the  natural  frequency  of  the  system  must  result,  there¬ 
fore,  from  a  pitch  moment  about  the  virtual  mass  center  caused  by  wave-in¬ 
duced  horizontal  forces  acting  on  the  four  legs  at  a  location  close  to  the 
undisturbed  free-surface  level. 

We  shall  here  assuem  that  the  pitch  springing  exciting  moment  for 
the  TLP  at  or  near  the  natural  frequency  can  be  obtained  from  the  springing 
excitation  for  a  single  leg.  In  fact,  we  shall  assume  that  the  pitch  ex¬ 
citing  moment  with  respect  to  the  virtual  mass  center  is 

F5  -  4(F*  *  zcFp  (3-1) 

where  F*  =  Fp^  +  f£^  =  pitch  springing  moment  for  a  single  leg  with 

respect  to  a  coordinate  system  located  at  the 
undisturbed  free  surface. 

F*  =  F^  +  f!|^  =  surge  springing  force  for  a  single  leg. 

zc  =  vertical  distance  from  undisturbed  free  sur¬ 
face  to  the  virtual  mass  center. 
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The  factor  "4"  in  equation  (3-1)  has  been  used  assuming  Lii.it  the  springing 
forces  for  the  four  legs  are  in  phase.  This  is  a  conservative  estimate. 
Figure  3-2  shows  schematically  the  springing  excitation  in  the  two  refer¬ 
ence  systems. 


Figure  3-2.  Total  TLP  pitch  moment  about  effective 
pitch  center. 

A  formal  derivation  of  the  first-  and  second-order  springing 
forces  are  presented  in  Salvesen,  et  al .  (1983).  It  is  shown  that  both 
the  first-  and  the  second-order  forces  can  be  expressed  in  terms  of  the 
first-order  velocity  potential,  4/^.  The  contribution  from  the  second- 
order  potnetial ,  ,  can  be  shown  to  be  small  for  the  vertical  cylinder  case, 

by  applying  an  approach  similar  to  the  method  used  by  MacCamy  and  Fuchs  (1954) 
and  by  assuming  that  the  waves  are  short,  A/H  «  1,  and  that  the  water  depth 
is  large,  A/d  <<  1,  it  can  be  shown  that  the  first-order  velocity  potential, 
for  the  vertical  cylinder  case  can  be  expressed  in  terns  of  a  series 
expansion.  Then  by  using  this  expansion  it  can  be  shown  that  the  first- 
a r. d  second  order  surge  forces  and  pitch  moments  can  be  expressed  in  terms 
of  :  '.o  Henkel1  s  function.  By  applying  this  approach  we  have  been  able  to 
develop  very  efficient  computational  methods  for  predicting  the  wave- 
inc  .Ci.d  springing  exciting  forces  and  moments. 
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Regular  Wave  Results  for  TLP  #1 


3.2.1  Pitch  Springing  Moment  Results 

In  this  section  we  shall  present  wave-induced  pi tch  spri ngi ng 


moment  results  for  the  TLP  #1  design 
The  first-order  springing  moment  is 

F^(t)  =  f^  cos  U)gt 

where  we  is  the  exciting  frequency, 
the  excitation  frequency  is  the  same 
order  springing  moment  is 

(t)  =  fg2^  cos  «et 

For  the  second-order  excitation,  the 
twice  the  wave  frequency.  Note  that 
mass  center  which  is  21.8  feet  below 


excited  by  regular  sinusoidal  waves. 

(3-2) 

Note  that  for  first-order  excitation, 
as  the  wave  frequency.  The  second- 

(3-3) 

excitation  frequency,  we>  is  equal  to 
the  pitch  moments  are  about  the  virtual 
the  undisturbed  free-surface  level. 


Figure  3-3  shows  both  the  first-crder  and  second-order  pitch 

springing  moment  amplitudes  as  function  of  the  exciting  period.  Results 

are  shown  for  two  wave  steepnesses,  H/A  =  0.050  and  0.075.  It  is  seen  from 

(2) 

the  figure  that  the  second-order  moment,  fl  ,  is  substantially  larger  than 

(1)  5 

the  first-order  moment,  f£  .  This  is  partly  due  to  the  fact  that  the  waves 
causing  the  first-order  moment  have  the  same  period  as  the  exciting  period, 

=  T^,  whereas  the  waves  causing  the  second-order  moments  have  periods 
twice  the  exciting  period,  Tw  =  2T^.  This  implies  that  the  wave  height  for 
the  waves  causing  the  second-order  moments  at  a  given  excitation  period  is 
tour  times  the  wave  height  for  the  waves  causing  the  first-order  moment 
(assuming  a  constant  H/A  value).  Furthermore,  it  is  seen  from  the  tabu¬ 
lated  values  in  Figure  3-3  that  at  the  natural  pitch  period  (T  =  2.506 
sec)  and  for  wave  steepness,  H/\  =  0.075,  the  first-order  moment  is  1600 
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3.2.2 


Equation  of  Motion 


If  we  consider  the  pitch  motion  about  the  virtual  mass  center, 
the  motion  is  uncoupled  and  the  equation  of  motion  is 

A'e  +  Be  +  ce  =  f5  •  e_1U)t  (3- 

Here  e  is  the  pitch  angle,  A  is  the  virtual  mass  moment  of  inertia,  B  is 
the  damping  coefficient  which  is  typically  a  function  of  the  pitch  angle, 

C  is  the  restoring  coefficient  and  fg  is  the  pitch  moment  amplitude.  As 
shown  in  Appendix  A,  for  the  TLP  #1  design 

A  =  4.06  *  1010  slug-ft2  (3-5 

C  =  2.553  *  1011  lb- ft  (3-6 

The  pitch  natural  period  is 

Tn  =  2tt  -^  =  2.506  seconds.  (3-7 

Of  interest  here  is  the  additional  tendon  tension  load  due  to 
pitch  springing  motion  rather  than  the  pitch  motion  itself.  The  tendon 
tension  load  per  tendon  is  linearly  proportional  to  the  pitch  angle 

t  =  aO  (3-8 

where  the  coefficient 

a  =  -ij  •  C33  ~  =  1.45  *  105  kips/rad.  (3-9 

Therefore,  the  tendon  tension  load  per  tendon  can  be  expressed  as 
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l(to)  = 


(3-10) 


(-oj2A+C)-iwB 


The  tendon  tension  load  can  also  be  expressed  in  terms  of  the  magnifica¬ 
tion  factor,  n,  as 


t(w)  =  n( fg 


(3-11) 


where  f^  a/C  is  the  zero  frequency  tendon  tension.  The  magnification  fac¬ 
tor  is  defined  as 

1 

n  =  - 

1  -  (-— )2  -  i2q  —  (3-12) 

“n  "n 

where  the  damping  factor  is  by  definition 


C 


B 

2fl“„ 


Here  Bc  is  the  critical  damping  coefficient. 


(3-13) 


Note  that  when  the  TIP  is  excited  by  a  pitch  moment  at  the  natural 
frequency,  u  ,  the  amplification  factor  is 


n 


1 

-i2C 


(3-14) 


3.2.3  Tendon  Springing  Loads 

By  equations  (3-11)  and  (3-14)  if  follows  that  the  additional 
tender,  tension  . oac  per  tendon  is 

tU  )  =  n  f,  £  =  -5  2.84  *  10"4  kips  (3-15) 

f.  a  L.  q 
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when  the  TLP  is  excited  at  the  natural  period,  =  2.506  seconds.  Here, 
the  pitch  exciting  moment  amplitude,  is  in  kips-ft  units.  Note  that  if 
the  damping  factor,  £,  is  a  constant  independent  of  the  pitch  motion  ampli¬ 
tude,  the  tendon  load,  x,  is  directly  proportional  to  the  pitch  moment,  f^. 
For  the  TLP  in  the  high  frequency  range,  the  damping  is  mainly  due  to 

•  viscous  forces, 

•  wave  generation,  and 

•  structural  damping. 

In  this  study,  we  have  included  damping  only  due  to  viscous  forces  and  wave 
generation.  It  is  expected  that  the  structural  damping  is  of  the  same  order 
of  magnitude  as  the  viscous  damping;  however,  it  is  difficult  to  obtain  ac¬ 
curate  estimates  for  the  structural  damping  without  a  complete  knowledge  of 
the  details  of  the  design  of  the  tendons. 

Figure  3-4  shows  the  damping  factor,  C,  as  a  function  of  tendon 
tension  load,  x  (not  that  x  is  linearly  proportional  to  the  pitch  angle, 
equation  (3-9)).  Also  shown  in  Figure  3-4  is  the  tendon  tension  load,  x, 
as  a  function  of  the  pitch  moment,  f ^ .  These  curves  have  been  obtained  by 
a  computer  program  which  models  the  viscous  forces  for  high  frequency  pitch 
motions  as  well  as  includes  the  wave  generation  damping.  Note  that  Figure 
3-4  is  only  applicable  to  the  TLP  §\  design  when  pitching  at  the  natural 
period,  Tn  =  2.506  seconds. 

Example  Results 

Consider  first-order  regular  wave  excitation  at  the  natural 
period  and  assume  that  H/A  =  0.075.  From  Figure  3-3  it  follows  that  the 
first-order  pitch  moment  amplitude  is 


Simplified  Prediction  Method 


The  contribution  to  the  variance  (i.e.,  the  square  of  the  RMS 
value,  a2)  for  second-order  excitation  can  be  estimated  by  the  following 
simple  procedure. 


w 

=/|HtU)|2  •  du> 

0  oo 

/  |Ht("'>2  d" 


.  (2)  /  \  /Ou2  TT  1 

=  (“nHA>  4  to 


m 


C 


(3-36) 


Applying  this  relationship  to  the  irregular  sea  condition  con¬ 
sidered  here,  we  have  approximately  the  same  result  as  obtained  by  numeri¬ 
cal  integration,  namely 

o2  ~  1980  . 

3.5  Sea  State  Investigation 

In  this  section  we  shall  present  irregular  sea  results  for  TLP  #1 
using  a  fully  developed  (FD)  Pi erson-Moskowi tz  sea  spectrum  as  well  as  the 
partially  developed  (PD)  Bretschneider  storm  sea  used  in  Section  3.4.  For 
both  cases  the  spectrum  is  given  by  equation  (3-34)  with  h^3  =  10  feet. 

For  the  FD  case 

T  =  8.73  seconds 
o 

and  for  the  PD  case  TQ  is  the  same  as  used  in  Section  3.4,  namely, 

T  =  7.04  seconds, 
o 

Figure  3-10  shows  the  sea  spectra  for  the  two  cases  considered, 
and  Figure  3-11  shows  the  second-order  pitch  exciting  moment  power  spectrum, 
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the  high  frequency  range. 


2 

Table  3-2  shows  the  square  of  the  RMS  values,  oT  (see  equation  3-28)  for  the 
first-  and  second-order  tendon  tension  response  spectrum  for  the  high-  and 
low-frequency  ranges  separately. 

The  results  presented  in  Table  3-2  seem  to  demonstrate  that  the 
first-order  contribution  in  the  low  frequency  range  dominates  the  total 
contribution  to  the  RMS  value  for  the  tendon  tension  response.  Considering 
only  the  low  frequency  range  contribution,  we  find  that  the  significant 
tendon  tension  value  for  this  sea  condition  is 

t  =  2.0  o  =  75  ki  ps  . 

If  we  consider  only  the  high  frequency  range,  we  have 
Xl/2  =  2.0  a  =  102  kips. 

It  should  be  stressed  that  the  main  objective  with  this  investi¬ 
gation  has  been  to  estimate  the  high  frequency  tendon  tension  load.  The 
low  frequency  computations  were  only  included  to  demonstrate  that  the  low 
frequency  tension  loads  may  be  as  important  as  the  high  frequency  part. 


Area  and  RMS  values  for  tendon  tension 
response  spectrum,  o2  and  o  in  kips2 
and  kips  1 


High  Frequency  Range 

Low  Frequency  Range 

Area 

RMS 

Area 

RMS 

Firs  —  Orrier  Contribution 

610 

24.7 

1420 

37.6 

Second-Order  Contribution 

1980 

44.5 

0 

0.0 

Total 

2590 

50.9 

1420 

37.6 

Table  3-2.  First-  and  second-order  contributions  to  area  and  RMS 
values  for  the  tendon  tension  response  spectrum  for 
the  high  and  low  frequency  ranges. 


c  =  0.05  *  10 

Figure  3-8  shows  the  tendon  tension  RAO,  Ht(w),  for  this  assumed  constant 
damping  factor. 
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Figure  3-8. 
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Tendon  tension  RAO,  HT,  for  TLP  #1, 
wn  =  2.508  rad/sec  and  damping  factor 
C  =  .0005 


Figures  3- 9a  and  3-9b  show  the  first-  and  second-order  tendon  ten 
sior,  response  power  spectrums.  As  indicated  :n  these  figures,  the  response 
power  spectrum  car  conveniently  be  divided  into  two  parts: 


the  low  frequency  range,  and 


3.4.2 


Computed  TLP  Results 


In  this  section  we  shall  present  computed  results  for  the  TLP  #1 
design  in  an  irregular  seaway  represented  by  a  Bretschneider  energy  spectrum 
for  a  partially  developed  storm  sea 

$  (w)  =  4  e"B/u)4  (3-34) 

a  o 

O) 

2  4  4 

where  A  =  486  hy^/^Q  ar>d  B  =  1948.18/Tq.  For  the  example  case  presenter 
here 


h.  ..  =  10  feet  and  T  =  7.04  sec. 
1/3  o 


Figure  3-7aand  3- 7b  show  the  first-  and  second-order  pitch  exciting  moment 
power  spectrum  as  obtained  by  equations  (3-23)  and  (3-25).  Please  note 
that  for  this  investigation,  the  first-order  pitch  moment  RAO  values, 


h[^(u)),  have  been  estimated  by  using  single  vertical  leg  results,  the 


approach  d’scussed  in  Section  3.1.2.  This  approach  is  only  valid  in  the 
high  frequency  range.  Therefore,  any  results  shown  here  in  the  low  frequency 
range  should  only  be  considered  as  an  order  of  magnitude  estimate.  It  is 


seen  in  Figure  3-7a  that  in  the  low  frequency  range  the  first-order  exciting 


moment  power  spectrum,  <+>'  ‘  is  as  expected,  an  order  of  magnitude  larger 


than  the  second-order  spectrum,  <|> 


(2) 


However,  as  seen  in  Figure  3- 7b  in 

(2) 


the  high  frequency  range,  near  the  pitch  natural  frequency,  (j>  '  is  sub¬ 


stantially  larger  than  $ 


(1) 


The  tender,  tension  load  response  power  spectrum  is  given  by  equa¬ 
tion  (3-26)  and  is 

‘^(w)  =  |HxU)  |2  •  (<J>^^(w)  *-  ^2/(w))  (3-35) 

where  is  the  tendon  tension  RAC  For  simplicity  we  shall  assume  that 

for  tne  'rr>_gu'ar  wave  results  presented  here  the  damping  factor  is 


and  the  average  of  the  1/10  highest  amplitudes  of  the  tendon  tension  load  is 


x  1/10 


2.55  0  . 

T 


(3-31) 


In  predicting  the  RMS  value  for  the  tendon  tension  load,  it  is 
important  to  recognize  that  the  pitch  damping  coefficient  is  nonlinear.  It 
is  here  suggested  that  the  nonlinear  damping  can  be  adequately  represented 
by  using  an  equivalent  linearization  method  where  the  damping  coefficient 
is  assumed  constant  and  equal  to  the  same  value  as  for  harmonic  motions 
at  the  natural  frequency  with  motion  amplitude  equal  to  the  significant 
motion  amplitude  for  the  irregular  wave  case. 


In  other  words,  the  following  procedure  is  used  in  predicting 
the  tendon  springing  load  in  an  irregular  seaway. 


First  assume  a  value  for  the  significant  tendon  tension  load, 

+ 

which  we  shall  label  T-^.  T*ien  using  this  value,  the 
first  estimated  damping  coefficient  is  obtained  from  Figure 
3-4.  The  significant  tendon  tension  load  is  then  computed 
following  the  above  steps  with  the  estimated  damping  coeffi¬ 
cient.  This  procedure  is  repeated  until  satisfactory  accu¬ 
racy  has  been  obtained. 


Note  that  if  the  free-surface  elevation  for  the  unidirectional 
seaway  is  given  by 

y(t)  =  2Xe'1Wnt  (3-32) 

n 

then  the  second-order  pitch  springing  moment  as  a  function  of  time  can  be 
expressed  as 

F52)U)  =  2  LAMV*(V  wm}  e'i(a>n  +  U>m)t  •  (3-33) 

n  m 
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The  pitch  moment  power  spectrum  due  to  second-order  wave- induced  excita¬ 
tion  is 

co/2 

=  2  •J' <J>a(|+  Aco)  cj)a  (^-A w) 

o  (3-25) 

•  ^-Aco)|^  d(Aoo) 

where  J  (u^,  002)  is  the  second-order  two-wave  pitch  springing  moment  per 
unit  wave  amp! itude. 


The  tendon  tension  load  response  power  spectrum  is 
4>t(w)  =  |Ht(o>)  j 2  •  (♦£«(«)  +  u)))  . 


(3-26) 


where  the  tendon  tension  per  unit  pitch  moment  excitation  is  defined  by 

(3-27) 


H  (U)  -  iM 
T  f5 


and  where  t(u>)  is  given  by  equation  ( 3- 10 ) .  The  variance  (i.e.,  the  square 
of  the  RMS  value)  of  the  tendon  tension  load  is 


-  I  <b  (u)  dw  . 
T  J  T 


(3-28) 


The  average  tendon  tension  load  amplitude  is 


i  =  1.25  •  a  , 
avg  t 


(3-29) 


:he  significant  value  of  the  tendon  tension  load  amplitude  i: 


r.  . 


.? .  00  •  a 


(3-30) 


From  Figure  3-5  it  is  seen  that  the  RAO  value  for  the  second-order  sum-fre- 
quency  moment,^  (w^,  u^)  =  260  kips-ft/ft  for  -  1.054  and  to^  =  1-454  rad/ 
sec.  Hence,  the  amplitude  of  the  second-order  springing  moment  component 
with  exciting  frequency  equal  to  the  natural  frequency  is 

2  [f£2)j  =  2A1A2J(u)1,  u2)  =  3620  kips-ft. 

From  Figure  3-4  it  follows  that  the  tendon  tension  load  for  the  above  two- 
wave  case  is 


f  21 

t'  '  =  780  kips  per  tendon. 

3.4  Irregular  Wave  Results  for  TLP  #1 

3.4.1  Theoretical  Formulation 

We  shall  assume  that  the  seaway  is  unidirectional  and  represented 
by  a  known  power  spectrum 

<J>,U)  (3-21) 

a 

and  that  the  wave-induced  pitch  exciting  moment  has  a  power  spectrum  given  by 

4>>)  =  41}(“)  +  ^2)  +  *  *  •  (3-22) 

Here  the  pitch-moment  power  spectrum  due  to  fi rst-order  wave-induced  exci¬ 
tation  is 

41}  =  IH^,1)(m)|  2  •  *aU)  (3-23) 


where  the  first-order  pitch  exciting  moment  amplitude  per  unit  wave  ampli¬ 
tude,  RAO,  is 


WAVE  FREQUENCY,  RAD/5EC 


Figure  3-6.  Two-wave,  second-order  pitch  moment  per 
unit  wave  amplitude,  ^  (u^,  u^). 

Two-Wave  Example  (see  Figure  3-5) 

Consider  two  waves  with  frequencies 

ujj  -  1-054  and  u ^  =  1.454  rad/sec. 

Tne  sum  frequency  of  these  two  waves  is  equal  to  the  natural  frequency  for 

the  TLP  -1 

ijj  -  m-i  +  coo  =  ?.503  rad/sec 
n  i  4 

\s s u.r.t  .  nat  they  both  have  the  same  wave  slope 


2.jf(2)J  e-i(«i  +  u2)t  =  2  .  AlA2J(Ul,  o>2)  e"i(a)l  +  u,2)t  (3-; 


where  and  are  the  wave  amplitudes  for  the  two  waves  with  frequencies 
u)j  and  respecti vely.  Here  to^)  is  defined  as  the  two-wave  pitch 

springing  moment  amplitude  per  unit  wave  amplitude  (RAO). 

Figure  3-5  shows  J  (u)j,  for  TLP  #1  for  sum  frequencies  equal 
to  the  pitch  natural  frequency: 

U1  +  w2  =  “n  =  2.508  rad/sec. 

Figure  3-6  shows  u^)  for  the  four  following  sum  frequencies: 

uj  +  =  2.0,  2.5,  3.0  and  3.5  rad/sec. 

Recall  that  for  this  component  of  the  second-order  pitch  moment,  the  exci¬ 
tation  period  is  equal  to  the  sum  frequency. 


EXCITING  FREQUENCY 

260.0 

U)r  =  U)l  ♦  U)o  =  U> 

WAVE  FREQUENCY.  RAD/SEC 

Figure  3-5.  Two-wave,  second-oder  pitch  moment  per 
unit  wave  amplitude  for  oji  +  002  =  con  = 
2.508  rad/sec. 


The  second-order  springing  exciting  force  is  given  by 


(3-17) 


Note  that  the  amplitudes  for  the  double  frequency  components,  f.  1  and 
I  fp)  1  1  J  1 li 

I  J  >22  dre  1"cientl'ca^  to  i tude  for  the  single  regular  wave  case. 

The  amplitude  of  the  sum  frequency  component  of  the  second-order  springing 

force  which  is  due  to  interactions  between  the  two  wave  components  is 


2M2)]  =2fJ.2)(u>1,  u>2)  (3-18) 

By  following  a  formal  procedure  as  used  for  the  single-wave  case,  it  can  be 
shown  that  the  amplitude  of  the  sum  frequency  component  (3-18)  can  be  ex¬ 
pressed  in  terms  of  the  first-order  potentials  for  the  two  individual  waves, 

1  and  ^  .  By  expressing  these  potentials  in  terms  of  series  expansions 

a  very  efficient  computer  program  has  been  developed  for  computing  the  two- 
(2) 

wave  amplitude,  f  1  J(wj,  ojg )  for  a  single  vertical  cylinder.  The  total  pitch 

moment  for  the  two-wave  case  is  obtained  by  equation  (3-1).  Note  that  a 

(21 

term  containing  the  second-order  potential,  »  has  been  ignored  in  this 
approach.  It  is  believed  that  this  term  will  result  in  a  very  small  contri¬ 
bution  to  the  total  force. 


Now  consider  two  regular  waves  with  sum  frequency  equal  to  the 
natural  frequency 

u)n  =  +  ui2  (3-19) 

tr.or  the  second-order  sum  frequency  component  of  the  pitch  moment  will  have 
••A.itation  frequency  identical  to  the  natural  frequency  and  may  be  expressed 
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If  we  assume  that  H/X  =  0.075,  it  follows  from  Figure  3-3  that  the  pitch 
moment  amplitude  is 

f[2^  =  5680  kips- ft. 

From  Figure  3-4  we  have  that  the  tendon  tension  load  is 
(2) 

t  '  =  1020  kips  per  tendon. 

The  above  results  for  first-  and  second-order  regular  wave  exci 
tation  are  summarized  in  Table  3-1. 


Pitch  moment 
amplitudes  in 
kips-ft 

Tendon  tension 
loads  in  kips 

First-order 
regular  wave 
excitation 

H/A  =  0.075 

1600 

470 

Second-order 
regular  wave 
excitation 

H/A  =  0.075 

5680 

1020 

Table  3-1.  Pitch  moment  amplitudes  and  tendon  tension 
loads  due  to  first-  and  second-order  regu¬ 
lar  wave  excitation. 

3.3  Two-Wave  Results  for  TLP  #1 


(2) 

(co)  for  these  two  sea  spectra.  It  is  seen  in  Figure  3-11  that  the  pitch 

exciting  power  spectrum  value  at  the  pitch  natural  frequency  (<a  =  2.508  rad/ 

sec)  is  almost  three  times  as  large  for  the  PD  sea  as  it  is  for  the  FD  sea 

condition.  According  to  equation  (3-36),  this  implies  that  the  variance, 

2 

o  ,  should  also  be  approximately  three  times  larger  for  the  PD  case  than  for 
the  FD  case. 

Table  3-3  shows  the  first-  and  second-order  contribution  to  the 
RMS  values  for  the  tendon  tension  responses  for  the  high-  and  low-frequency 
ranges  both  for  the  RD  and  the  PD  sea  states.  It  is  seen  that  for  the 
high  frequency  range  the  first-  and  second-order  RMS  values  for  the  PD  sea 
state  are  about  v^times  as  large  as  they  are  for  the  RD  sea  state.  Whereas 
for  the  low  frequency  range  where  only  the  first-order  contribution  is  of  a 
significant  magnitude,  the  RMS  values  for  the  FD  sea  state  are  larger  (by  a 
factor  of  two)  than  those  for  the  PD  sea  state. 

These  results  seem  to  indicate  that  it  is  important  to  consider 
both  partially  developed  and  fully  developed  sea  states  when  predicting  the 
tendon  tension  responses  due  to  wave-induced  springing.  In  general,  it  can 
be  expected  that  the  second-order  contribution  will  be  substantially  larger 
in  a  partially  developed  sea  state  than  in  a  fully  developed  sea  state  with 
the  same  significant  wave  height. 


RMS  values  for  the  tendon  tension 
response  spectrum,  in  kips 

High  Frequency  Range 

Low  Frequency  Range 

FD  Sea 

PD  Sea 

FD  Sea 

PD  Sea 

First. -Order  Contribution 

16.1 

24.7 

75.7 

37.6 

Second-Order  Contribution 

26.3 

44.5 

0.1 

0.0 

Total 

3C.9 

50.9 

75.7 

37.6 

Table  3-3.  First-  and  second-order  contributions  to  RMS  values  for 
the  tendon  tension  response  spectrum  for  the  high  and 
low  frequency  ranges  for  fully  developed  ( FD)  and  par¬ 
tially  developed  (PD)  sea  states. 


3-26 


In  this  section  we  shall  investigate  the  effect  of  the  location 
of  the  pitch  center  for  TLP  #1.  Note  that  the  pitch  center  is  defined  as 
the  virtual  mass  center  (including  both  the  mass  of  the  TLP  and  the  hydro- 
dynamic  added  mass).  We  shall  consider  two  pitch  center  locations:  the 
actual  location  for  the  TLP  #1,  namely 

z ^  =  -21.8  feet 
and  25"  of  this  distance, 

=  0.25  =  -5.5  feet. 

3.6.1  Regular  Wave  Results 

Figure  3-12  shows  the  first-  and  second-order  pitch  springing 
moments  for  TLP  HI  with  the  two  stated  pitch  center  locations.  It  is  seen 
from  this  figure  that  at  the  natural  period  the  second-order  pitch  moment  is 
reduced  by  a  factor  of  about  10  by  decreasing  the  pitch  center  to  25%  of 
the  design  value.  The  first-order  moment  is  reduced  by  a  factor  of  about  50 

Table  3-4  shows  the  oitch  moment  amplitude  values  at  the  natural 
period  as  obtained  from  Figure  3-12.  Also  shown  in  Table  3-4  are  the  tendon 
tension  loads  due  to  these  pitch  moment  values.  The  loads  have  been  obtaine 
using  the  nonlinear  relationship  between  the  pitch  moment  and  tension  loads 
as  given  by  Figure  3-4.  It  is  seen  from  the  tendon  tension  values  given  in 
Table  3-4  that  the  second-order  tendon  tension  load  is  reduced  from  1020 
kips  to  240  kips  by  decreasing  the  pitch  center  to  25%  of  the  design  value. 
Similarly,  the  first-order  tendon  tension  load  is  reduced  from  470  kips  to 


1600 

470 

30 

— 

24 

Second-Order  Single-Wave 
Excitation,  H/X  =  0.075 


Table  3-4.  Pitch  moment  amplitudes,  fc  and  tendon  tension  loads, 
t,  due  to  first-  and  second-order  excitation  at  the 
natural  period  for  vertical  pitch  center  at  zi  =  -2l.8  ft 
and  at  12  -  -5.5  ft. 


3.6.2 


Irregular  Wave  Results 


Figure  3-13  shows  the  power  spectra  for  the  second-order  pitch 
springing  exciting  moments  for  the  two  pitch  center  locations.  It  is  seen 
that  the  second-order  power  spectrum  has  been  reduced  by  an  order  of  magni¬ 
tude  by  decreasing  the  pitch  center  to  25%  of  the  design  value. 

Table  3-5  shows  the  RMS  values  for  the  tendon  tension  response  for 
the  two  pitch-center  locations.  It  is  seen  that  there  is  a  large  reduction 
in  the  tendon  tension  RMS  values  for  the  contribution  from  the  high  frequency 
range.  However,  for  the  low  frequency  range,  where  the  loads  are  dominated 
by  first-order  excitation,  there  is  actually  a  substantial  increase  in  the 
tendon  tension  when  the  pitch  center  is  decreased  to  25%  of  the  design  value. 
The  reason  for  this  difference  between  the  high  frequency  and  the  low  fre¬ 
quency  range  is  that  the  high  frequency  range  excitation  is  caused  by  short 
waves  with  exciting  force  locations  close  to  the  free  surface,  whereas  the 
low  frequency  excitation  is  caused  by  longer  waves  with  exciting  force  loca¬ 
tions  below  the  design  pitch  center.  Hence,  in  the  low  frequency  range 
the  pitch  exciting  moments  are  actually  increased  when  the  pitch  center  is 
moved  closer  to  the  free-surface. 

In  summary,  the  pitch  center  investigation  results  show  that  a 
very  effective  way  to  reduce  the  second-order  tendon  tension  springing  loads 
is  to  move  the  pitch  center  closer  to  the  free  surface  level.  However,  the 
first-order  tendon  loads  in  the  low  frequency  range  will  most  likely  increase 
by  moving  the  pitch  center  closer  to  the  free  surface.  Therefore,  it  is 
very  difficult  to  determine  if  it  is  advantageous  to  move  the  pitch  center 
location.  A  careful  investigation  of  the  first-  and  second-order  tendon  load 
results  for  both  the  high  and  the  low  frequency  ranges  is  required  before  any 
judgment  with  regard  to  the  pitch  center  location  can  be  made. 
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Figure  3-13.  Second-order  pitch  exciting  moment  power  spectra 
for  pitch  center  zc  =  -21.78  ft  and  iQ  -  -5.45  ft. 


RMS  values  for  the  tendon  tension 
response  spectrum,  in  kips 

equency  Range  Low  Frequency  Range 


High  Frequency  Range 


-5.5 


z\  =  -21.8'  z. 


First-Order  Contribution 


Second-Order  Contribution 


Total 


Table  3-5.  First-  and  second-order  contributions  to  RMS  values  for  the 
tendon  tension  response  spectrum  for  the  high  and  low  fre¬ 
quency  ranges  for  vertical  pitch  center  at  zi  =  -21.8  ft  and 
Zj =  -5.5  ft. 


Zl  =  -21.8* 

z  £  =  -5.5' 

37.6 

59.6 

0.0 

0.0 

37.6 

59.6 

In  this  section  we  shall  compare  springing  results  for  the  TLP  #1 
and  TLP  ft 2  designs.  A  detailed  description  of  the  two  TLPs  is  included  in 
Appendix  A.  Note  that  TLP  ft 2  is  1.2  times  larger  than  TLP  ffl  and  that  TLP 
HI  is  assumed  to  be  designed  for  20Q0- foot  depth  whereas  TLP  # 2  is  designed 
for  3000- foot  depth. 

3.7.1  Regular  Wave  Results 

Figure  3-14  shows  the  first-  and  second-order  pitch  moment  for  TLP 
ffl  and  TLP  ft?  in  regular  waves  with  wave  steepness,  H/A  =  0.075.  It  is 
seen  that  the  pitch  moments  are  somewhat  larger  for  the  larger  TLP.  However, 
more  importantly  is  the  fact  that  the  natural  period  for  the  larger,  deep¬ 
water  TLP  #2  is  3.134  sec,  whereas  for  the  smaller  TLP  ffl  it  is  only  2.506 
sec.  The  longer  waves  exciting  TLP  #2  at  its  natural  period  have  much  more 
energy  than  the  shorter  waves  exciting  TLP  #1  at  its  natural  period. 

The  numerical  values  for  the  pitch  moment  amplitudes  and  the 
resulting  tension  loads  at  the  natural  periods  are  presented  in  Table  3-6. 

It  is  seer,  that  the  second-order  pitch  moment  for  TLP  #2  is  about  three 
times  as  large  as  for  TLP  ffl  and  that  the  resulting  tendon  loads  for  TLP  #2 
are  about  twice  as  large  as  for  TLP  #1. 

These  results  seem  to  indicate  that  the  springing  problem  is  much 
more  critical  for  a  deep-water  TLP  with  a  higher  natural  period. 


3.7.2 


Irregular  Wave  Results 


Figure  3-15  shows  the  power  spectrum  for  the  second-order  pitch 
springing  exciting  moment  for  both  TLP  #1  and  TIP  #2.  It  is  seen  that  at 
the  natural  frequencies  for  the  two  platforms  the  value  of  the  power  spec¬ 
trum  is  about  twice  as  large  for  TLP  #2  as  for  TLP  #1. 

Table  3-7  shows  the  RMS  values  for  the  tendon  tension  response  for 
both  TLP  HI  and  TLP  i 2 .  It  is  seen  that  the  RMS  values  are  somewhat  larger 
for  the  larger  TLP,  but  not  as  much  larger  as  the  springing  exciting  moments 
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Figure  3-15.  Second-order  pitch  exciting  moment  powe1''  spe 
for  TLP  #1  and  TLP  #2. 


RMS  values  for  the  tendon  tension 
response  spectrum,  o  in  kips 

High  Frequency  Range 

Low  Frequency  Range 

TLP  #1 

TLP  #2 

TLP  #1 

TLP  #2 

First-Order  Contribution 

24.7 

71.9 

37.6 

40.5 

Second-Order  Contribution 

44.5 

51.8 

0.0 

0.0 

Total 

50.9 

88.6 

37.6 

40.5 

Table  3-7.  First-  and  second-order  contributions  to  RMS  values  for 


the  tendon  tension  response  spectrum  for  the  high  and  low 
frequency  ranges  for  TLP  #1  (X  =  1.0)  and  TLP  $2  (X  =  1.2' 
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REDUCED  NUMBER  OF  TENDONS  INVESTIGATION 


4.1  Objective 

The  main  objective  of  this  part  of  the  project  has  been  to  inves¬ 
tigate  the  surge  motions  and  the  tendon  tension  loads  for  the  case  of  tendon 
failure,  i.e.  only  three  tendons  for  one  of  the  legs. 


4. 2  Non-Rigid  TLP  Structure 

First,  it  shall  he  assuemd  that  the  TLP  is  not  completely  rigid 
so  that  aJ_T_  of  the  loads  from  the  broken  tendon  are  taken  up  by  the  three 
remaining  tendons  at  that  leg. 

4.2.1  Maximum  External  Loading 


If  we  let  tQ  equal  the  pretension  value  per  tendon  for  a  normal 
operational  condition,  then  the  pretension  value  for  each  of  the  three 
remaining  tendons  is 


^  i-33  V 


(4-1) 


We  shall  use  an  asterisk  to  indicate  that  the  value  is  for  the  case  with 
one  tendon  removed.  We  shall  assume  that  the  tendons  are  designed  to  take 
an  additional  loading  equal  to  the  pretension  value.  This  implies  that  the 
maximum  total  loading  at  the  leg  with  three  tendons  is 


3t. 


max 


V 


(4-2) 


Subtracting  the  pretension  value  f  r  these  three  tendons  we  have  that  the 
total  additional  loading  which  can  be  taken  up  by  these  three  tendons  is 


6to 


(4-3) 
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Since  the  total  additional  maximum  loading  at  a  leg  with  four  intact  legs  is 


4t 


max 


(4-4) 


it  is  seen  by  comparing  (4-3)  and  (4-4)  that  the  leg  with  three  tendons  can 
only  experience  half  the  external  loads  as  a  leg  with  four  tendons  if  it 
is  required  that  the  total  load  per  tendon  shall  not  exceed  a  value  which 
is  twice  the  pretension  value. 


4.2.2  Surge  Motion  Characteristics 


If  we  consider  the  coefficients  and  the  exciting  forces  for  the 
surge  equation  of  motion  (see  equation  (3-1))  it  is  easily  recognized  that 
the  mass,  added  mass  and  damping  coefficient  will  remain  about  the  same  with 
one  less  tendon.  Similarly,  the  excitation  forces  will  be  about  the  same. 
There  will  be  some  reduction  in  the  viscous  drag  due  to  the  elimination  of 
one  tendon;  however,  this  effect  is  small  compared  to  the  total  viscous 
force.  Also  note  that  the  (linear)  surge  restoring  coeffiicent  which  is 
given  by 


C 


11 


a 


"t 

T 


(4-5) 


will  be  approximately  the  same.  Here  T0,  which  is  the  total  pretension 

value,  is  unchanged  and  w^,  which  is  the  weight  of  the  tendons,  is  only 

reduced  by  about  6%.  Note  also  that  the  wt/2  term  is  small  relative  to 

the  ~  term . 
o 

It  can  be  concluded,  therefore,  that  the  surge  motion  characteris¬ 
tics  of  the  TIP  will  remain  practically  the  same  with  one  less  tendon.  No 
additional  computations  are  needed  to  confirm  this. 
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4.2.3 


Springing  Loads 


If  we  consider  the  heave  restoring  coefficient  separately  for  each 
of  the  four  legs,  we  have  that  the  heave  restoring  coefficient  for  the  leg 
with  three  tendons  is  75%  of  the  restoring  coefficient  for  the  other  legs, 
hence 


*C3®g  =  0.75  C339.  (4-6) 

This  implies  that  the  natural  heave  period  for  the  leg  with  three  tendons  is 

T*  =  Tn/y05  =  1.15  •  Tn  (4-7) 

where  Tn  is  the  natural  heave  period  for  the  other  three  legs.  Hence,  the 
leg  with  only  three  tendons  has  a  15%  larger  natural  period  than  the  other 
legs. 


The  springing  laods  are  caused  by  vertical  motions  at  the  individ¬ 
ual  legs  resulting  from  natural  frequency  pitch  motions  of  the  total  TLP. 
With  three  legs  having  a  natural  period  of  Tn  and  one  leg  having  a  natural 
period  equal  to  1.15  T  ,  it  is  clear  that  the  pitch  motion  will  not  have  a 
single  distinct  natural  pitch  period.  Therefore,  it  can  be  concluded  that 
the  springing  load  problem  which  is  mainly  due  to  a  very  distinct  natural 
frequency  will  not  be  as  critical  as  it  is  for  the  case  with  all  of  the 
tendons  intact. 

4.3  Rigid  TLP  Structure 

It  shall  now  be  assumed  that  the  TLP  structure  is  completely 
rigid  and  that  the  load  from  the  broken  tendon  is  taken  up  by  all  of  the 
remaining  15  tendons.  If  we  use  a  numbering  system  as  shown  in  Figure  4-1 
where  leg  #1  is  the  leg  with  only  three  tendons,  we  have  that  the  pretension 
value  is  increased  as  shown  below: 
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Three  tendons  at  leg  #1 
Eight  tendons  at  legs  #2  and  #3 
Four  tendons  at  leg  #4 


Pretension  Increase 
25% 


6.25% 

-6.25% 

In  this  case  the  maximum  external  load  which  the  total  platform 
can  experience  is  56%  of  the  load  if  can  experience  with  all  of  the  tendons 
intact  if  it  is  assumed  that  the  maximum  design  load  per  tendon  is  twice 
the  pretension  value.  This  result  was  obtained  by  recognizing  that  the 
three  tendons  at  leg  #1  can  only  experience  an  additional  load  given  by 

3(0.75  tQ)  =  2.25  tQ  (4-8) 

whereas  four  intact  tendons  are  designed  for  4  t  additional  load. 

Again,  for  the  rigid  TLP  case  the  surge  motion  characteristics 
will  be  practically  unchanged  by  eliminating  one  single  tendon. 

The  springing  problem  in  this  case  becomes  a  little  more  compli¬ 
cated.  First,  it  can  be  shown  that  the  pitch  center  will  move  a  distance 
of  0.047  L  (where  L  is  the  leg  spacing)  as  shown  in  Figure  4-1.  Further¬ 
more,  it  can  be  shown  that  the  pitch  restoring  coefficient  for  the  case  of 
one  less  tendon  is 


C*5  =  0:864  C55  (4-9) 

where  is  the  restoring  coefficient  with  the  tendons  intact.  The  pitch 
natural  period  will  be 


1.076  Tn 


where  Tn  is  the  intact  pitch  period. 
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Consider  springing  tendon  tension  loads  caused  by  natural  fre¬ 
quency  pitch  motions  excited  by  second-order  wave- induced  forces.  Since 
the  second-order  excitation  is  proportional  to  the  wave  height  squared,  and 
the  wave  height  is  proportional  to  the  period  squared  (assuming  a  fixed  wave 
slope,  i.e.,  H/  equal  to  a  constant),  it  follows  that  the  second-order 
pitch  springing  moment  at  a  natural  period  equal  to  1.076  Tn  is  approxi- 
ma tely 

f*  =  (1.076)4  •  f5  =  1.34  f5  (4-11) 

where  f5  is  the  pitch  moment  at  a  natural  period  equal  to  T^.  In  other 
words,  the  pitch  springing  moment  in  this  case  is  about  34%  larger. 

For  natural  frequency  motions  with  all  of  the  tendons  intact,  the 
pitch  angle  is  given  by  (see  Section  3.2.2) 

6  =  f5/un  B.  (4-12) 

For  the  case  with  one  tendon  removed,  the  pitch  angle  is 

0*  =  f*/u*  B.  (4-13) 

b  n 

It  can  be  assumed  that  the  damping,  B,  remains  the  same.  Since  the  new 
natural  frequency  is 

a )*  =  0.929  un  (4-14) 

it  follows  from  equation  (4-13)  by  introducing  equations  (4-11)  and  (4-14) 
that 

9*  =  (1.34  •  f •  )/(0. 929  a*  )B  =  1.449. 

Hence  the  pitch  angle  will  be  44%  larger  than  for  the  intact  case.  The 
loading  in  the  tendons  is  proportional  to  the  pitch  angle 


where  the  proportionality  factor,  a,  for  the  case  with  one  less  tendon  is 
a*  =  1.067  •  j  a  =  1.42a. 

Therefore,  the  springing  tendon  tension  per  tendon  at  leg  #1  will  be 
T*  =  a*  •  0*  =  1.42a  •  1.440  =  2.04t. 

This  implies  that  if  we  consider  regular  wave  springing  excitation,  the  TLP 
with  one  less  tendon  will  experience  tendon  laods  at  the  leg  with  only  three 
tendons  which  are  2.04  times  the  springing  load  in  the  tendons  for  a  TLP 
with  all  of  the  tendons  intact. 

Now  considering  that  the  pretension  value  for  the  three  tendons 
at  leg  #1  has  already  been  increased  by  25%,  it  seems  as  if  the  additional 
loading  due  to  springing  can  become  quite  critical  for  the  case  with  one 
less  tendon. 
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Appendix  A 


PRINCIPAL  DIMENSIONS  AND  NATURAL  PERIODS  FOR 
TLP  #1  AND  TLP  #2 


Two  TLP  designs  are  considered  in  this  investigation.  They  are 
referred  to  as  TLP  #1  and  TLP  #2  with  the  following  principal  dimensions: 


tlp  n 

TLP  #2 

Vertical  Leg  Diameter,  D^ 

57.0  ft 

69.0  ft 

Pontoon  Diameter, 

31.0  ft 

37.0  ft 

Leg  Centerline  Spacing,  L 

220.0  ft 

260.0  ft 

Draft  of  Leg,  H 

121.0  ft 

145.0  ft 

Depth  to  Ocean  Floor,  d 

2000.0  ft 

3000.0  ft 

Note  that  the  two  platforms  have  the  same  geometry  with  a  scaling  factor 
of  approximately  1.2;  however,  the  depth  for  TLP  #2  is  1.5  times  the  depth 
for  TLP  #1.  The  two  TLPs  have  a  configuration  as  shown  in  Figure  A-l  with 
four  vertical  legs  connected  with  four  horizontal  pontoons  and  some  kind 
of  deck  structure.  In  addition  to  this  main  structure,  there  is  also  an 
appropriate  number  of  braces  for  which  volume  and  mass  are  included  by 
adding  an  estimated  percentage  factor.  In  the  hybrid-finite-element 
method  (HFEM)  computations  only  the  legs  and  the  pontoons  have  been 
included.  There  are  four  vertical  tendons  at  each  of  the  four  legs. 


« -  l - * 

7 - 1 - 1 - j- 


4  tendons 
at  each  leg 


Figure  A-l.  TLP  Geometry. 
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The  buoyancy  of  TLP  # 1  is 


Ag  =  116,880  kips 

and  the  total  pretension  value  (including  tendons  and  risers)  is 
Tq  =  40,100  kips 
with  the  weight  given  by 

WQ  =  Ag  -  Tq  =  76,700  kips. 

The  total  weight  of  the  tendons  and  risers  in  air  is 
w a  =  13,900  kips 

a 

and  in  water 

ww  =  12,700  kips. 

The  effective  mass  of  the  total  TLP  system  is 

w 

M  =  (Wo  +  f)  ?  =  2,60  *  1q6  slugs 

and  the  surge  restoring  coefficient  is 

w 

cll  =  {To  '  T]  1  =  18,32  kips/ft 

where 

H  =  d-  H-  e  =  1840  ft. 


Here  e  =  39  is  the  height  of  the  bottom  tendon  mounting. 


below  the  water  level  and  that  the  virtual  pitch  mass  moment  of  inertia 
about  this  center  is 


M55  +  A55  =  4.06  *  10-  slugs-fr 


and  hence  the  pitch  natural  period  is 

/M„  +  A,, 

T55  =  2irV  ~ c~  ' ~  =  2-506  sec- 

All  of  the  above  results  are  for  TLP  #1.  Table  A-l  shows  the 
results  for  both  TLP  #1  and  TLP  #2. 

Table  A-l.  Mass,  restoring  coefficients,  and  natural  periods  for 
HP  #1  and  TLP  #2. 


TLP  #1 

TLP  #2 

Buoyancy,  (kips) 

116,800 

202,000 

Pretension,  Tq  (kips) 

40,100 

69,000 

Weight,  W  (kips) 

76,700 

113,000 

Tendon  weight  in  air,  wfl  (kips) 

13,900 

23,200 

Tendon  weight  in  water,  w  (kips) 

W 

Mass,  M  (slugs) 

12,700 

21,900 

2.60  *  106 

4.44  *  106 

Surge  restoring  coef.,  (kips/ft) 

18.32 

20.7 

Surge  added  mass,  A^  (slugs) 

2.63  *  106 

4.28  *  106 

Surge  natural  period,  Tjj  (sec) 

106.12 

128.96 

Heave  restoring  coef.,  C33  (lb/ft) 

21.1  *  106 

24.05  *  106 

Heave  added  mass,  A33  (slugs) 

1.44  *  106 

2.64  *  106 

Heave  natural  period,  T33  (sec) 

Pitch  restoring  coef.,  (ft-lb) 

2.738 

3.37 

25.53  *  1010 

40.6  *  1010 

Pitch  center  location,  z„  (ft) 

21.8 

26.2 

c  0 

Pitch  virtual  mass,  M^^+Agg  (slg-ft^) 

4.06  *  1010 

10.1  *  1010 

Pitch  natural  period,  T^g  (sec) 

2.506 

3.13 
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DETERMINATION  OF  SWAY  DISPLACEMENT  COMPONENTS  BY 
SAI  AND  QUADRATIC  DRAG  FORMULAS 


B. 1  Computer  Run  Descriptions 

A  systematic  series  of  computer  runs  was  performed  with  the  non¬ 
linear  time-domain  surge  motion  computer  code.  The  individual  surge  dis¬ 
placement  components  were  determined  in  order  to  evaluate  the  differences 
in  each  component  using  the  SAI  and  quadratic  drag  formulas,  and  for  TLP  #1 
and  TLP  #2.  The  following  current,  wind  and  wave  conditions  were  selected 
for  this  investigation: 

(i)  Current  -  The  current  was  input  as  a  function  of  depth  below  the 
free  surface  as  follows: 


Z,  ft 

Current,  knots 

0. 

2.70 

121. 

1.70 

1000. 

0.60 

2000. 

0.30 

(ii)  Wind  -  Steady  wind:  W$  =  70  knots 

Unsteady  wind:  Wys  =  7  knots  cosoont 
where  wn  is  the  natural  frequency  of  the  system. 


(iii)  Waves  -  Maximum  design  wave:  H  =  80  ft,  T  =  14.0  sec 

For  this  wave,  the  primary  wave  exciting  force  and  phase  are 

F^  =  22,870  kips  and  =  1.596 
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The  steady  second-order  drift  force  is 


F 


(2) 


=  36  kips. 


For  the  given  wind  conditions,  with  steady  and  unsteady  components 
present,  the  wind  force  is 


Fw  =  960  kips  +  191  kips  cosa>nt. 

The  computer  runs  were  conducted  both  with  the  SAI  drag  formula 
(referenced  by  the  letter  "A")  and  the  quadratic  drag  formula  (referenced 
by  the  letter  "B").  The  following  runs  were  performed: 

1A  and  B  Current  alone  (no  wind  or  waves) 

2A  Maximum  design  wave,  primary  excitation  only  (no  drift  force, 

current  or  wind) 


3A  Maximum  design  wave  (primary  excitation  only),  with  inte¬ 

gration  up  to  the  disturbed  free  surfece  (no  drift  force, 
current  or  wind) 


4A  and  B  Current  and  maximum  design  wave,  primary  excitation  only 

(no  wind  or  drift  force) 


5A  and  B  Current,  maximum  design  wave  including  steady  drift  and 

steady  wind  (no  unsteady  wind) 

6A  Unsteady  wind  only  (no  current,  wave  or  steady  wind) 

7A  and  B  Current  and  unsteady  wind  (no  wave  or  steady  wind) 
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8A  and  B  Maximum  design  wave  (primary  excitation  only)  and  unsteady 

wind  (no  drift  force  or  steady  wind) 

9A  and  B  Current,  maximum  design  wave  including  steady  drift,  steady 

and  unsteady  wind 

The  conditions  stated  above  are  tabulated  in  Table  B-l. 


Table  B-l.  Conditions  used  for  surge-displacement  component 
investigation. 


*  Integration  up  to  free  surface 
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Discussion  of  Results 


The  numerical  results  for  the  cases  stated  above  are  given  in 
Table  B-2,  in  terms  of  the  steady  state  amplitudes  of  surge  at  the  wave 
frequency,  the  natural  frequency  and  the  steady  surge.  The  corresponding 
time-domain  plots  of  the  surge  motion  are  shown  in  Figures  B-l  through  B-9. 

Runs  1A  and  IB 

Figure  B-l  shows  the  surge  response  for  runs  1A  and  IB  with  cur¬ 
rent  only.  As  should  be  expected,  the  final  steady  displacement  is  the 
same  for  both  the  SAI  and  quadratic  drag  formulas, 

Xn  =  14.4  feet, 
o 

The  only  difference  between  the  results  for  the  two  current  only  cases  is 
the  transient  motions  die  much  slower  with  the  SAI  drag  than  the  quadratic. 
It  is  believed  that  the  quadratic  drag  dies  out  much  too  rapidly. 

Run  _2A 

Figure  B-2  shows  the  results  for  Run  2A,  regular  maximum  wave 
only  (no  drift  force).  It  is  seen  that  the  primary  motion  amplitude  is 

Xj  =  22.1  feet. 

Note  that  this  primary  motion  amplitude  remains  the  same  throughout  these 
computations  which  shows  (as  we  had  expected)  that  the  primary  motions  are 
not  in  any  way  effected  by  the  current  or  by  the  slowly-varying  motions. 
Furthermore,  viscous  damping  has  an  extremely  small  effect  on  these  motions 
(les$  than  1%)  and  does  not  need  to  be  included  in  primary  motion  predic¬ 
tions.  The  primary  motions  are  dominated  by  inertia  effects  and  hence  are 
linear  with  respect  to  wave  amplitude.  This  implies  that  linear  superposi¬ 
tion  and  linear  statistics  can  be  applied  to  the  primary  motions. 


Run  3A 


Figure  B-3  shows  the  results  for  Run  3A.  The  only  excitation  is 
the  primary  wave  excitation  (no  drift  force),  but  in  this  case  the  hydro- 
dynamic  pressure  is  integrated  over  the  instantaneous  wetted  surface  area. 

It  is  seen  that  in  this  case  there  is  an  additional  steady  sway  displacement 

XQ  =  7.7  feet. 

This  is  the  maximum  additional  displacement  due  to  integration  over  the 
instantaneous  wetted  surface  rather  than  over  the  undisturbed  wetted  sur¬ 
face.  The  difference  is  so  small  that  we  shall  not  include  this  effect  in 
any  of  the  other  computations. 

Runs  4A  and  4B 

Figure  B-4  shows  the  results  for  Runs  4A  and  4B,  current  and  pri¬ 
mary  wave  excitation.  It  is  seen  that  the  quadratic  damping  formula  predicts 
that  the  steady  displacement  is 


XQ  =  42.3  feet 

whereas  the  SAI  drag  formula  predicts  the  steady  displacement  to  be 

X  =16.1  feet, 
o 

The  difference  is  due  to  the  quadratic  damping  predicting  a  much  larger 
additional  drag  force  due  to  the  interaction  between  the  current  and  the 
primary  and  the  wave-particle  velocity  than  the  SAI  formula  prediction. 


In  summary,  the  quadratic  damping  predicts  that  the  displacement 
due  to  the  current  will  be  magnified  by  the  primary  motions  by  the  factor 


42^3 

14.4 


=  2.94 
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comparisons  for  the  linear  results.  It  should  also  be  kept  in  mind  tnat 
the  HFEM  computations  cost  at  least  $6,000  per  TLP  configuration  of  whicn 
about  $4,000  is  for  the  computer  time  and  $2,000  is  for  preparing  the 

input  data. 

For  the  "Nonlinear  Drift  Forces  and  Surge  Responses"  (Part  2) 
extensive  verification  only  exist  for  the  second-order  steady  drift  forces 
in  regular  waves  (Part  2a).  See  for  example  Figure  2  of  Salvesen  et  al . 
(1982)  where  three  different  theoretical  approaches  are  compared  with 
experimental  results.  Extensive  comparisons  between  the  drift  forces 
computed  for  TLP  configuration  and  experimental  results  have  been  pre¬ 
ferred  by  Chevren,  but  these  results  are  unpublished. 

There  is  a  great  need  for  verification  of: 

•  the  nonlinear  surge  response  (Part  2a),  and 

•  the  nonlinear  springing  forces  and  responses  (Part  3a  and  b). 

This  is  where  we  had  hoped  to  perform  some  valuable  new  comparisons  between 
theory  and  experiments.  But  there  is  almost  a  complete  lack  of  published 
experimental  data  for  controlled  environmental  conditions  for  which  computa¬ 
tions  can  be  performed.  In  order  that  a  useful  validation  can  be  performed, 
it  is  essential  that  not  only  precise  geometry  data  be  available,  but  also 
tnat  the  environmental  condition  be  accurately  known,  it  should  be  pointed 
out  that  we  have  not  been  able  to  locate  any  springing  data  whatsoever. 
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DISCUSSION  OF  COMPARISONS  BETWEEN 
THEORY  AND  EXPERIMENTS 

The  computed  TLP  results  presented  in  this  report  can  be  divided 
into  three  main  parts: 

Part  1:  Linear  Forces  and  Responses 

Linear  f requency-domai n  potential  flow  results  obtained  by  the 
HFEM  computer  program.  These  results  include  added-mass  coefficient, 
damping  coefficient  due  to  wave  generation  and  first-order  wave  exciting 
forces  and  moments  as  well  as  the  linear  wave-frequency  components  of  the 
motions. 

Part  2:  Nonlinear  Drift  Forces  and  Surge  Responses 

a)  Second-order  steady  draft  forces  for  regular  waves 

b)  Nonlinear  surge  responses. 

Part  3:  Nonlinear  Springing  Forces  and  Responses 

a)  Second-order  high-frequency  springing  forces  in  regular 
waves  and  in  wave  groups  consisting  of  two  regular  waves. 

b)  High  frequency  springing  response  load  in  the  tendexes. 

A  careful  verification  of  all  of  these  computed  results  would  be 
desirable  and  in  particular  it  would  be  useful  to  compare  such  results  to 
experimental  results.  Fortunately,  quite  satisfactory  comparisons  have 
already  been  completed  for  the  “Linear  Forces  and  Responses"  (Part  1). 

See  for  example  the  User's  Manual  for  the  HFEM  Computer  Program  (Yue,  1982). 
In  our  opinion  it  would  not  be  too  valuable  to  perform  extensive  additional 
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For  the  SAI  drag  formula  the  following  sectional  viscous  drag 
coefficients  were  used  for  all  of  the  computations  presented  in  this 
report.  (Note  that  the  drag  coefficients  were  reduced  due  to  sheltering 
effects  in  the  same  way  as  for  the  quadratic  drag  fomrula  discussed  on 
the  previous  page.) 


Member 

CDs 

r 

LD1 

CD2 

legs 

0.5 

0.5 

0.5 

Pontoons 

0.45 

0.5 

0.5 

Horizontal  Brace 

0.4 

0.7 

0.65 

Vertical  Member 

0.42 

0.7 

0.65 

Diagonal  Brace 

0.4 

0.7 

0.6 

Tendons:  Top 

0.7 

0.55 

0.5 

Botton  (1000  ft) 

1.2 

0.55 

0.5 
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VISCOUS  DRAG  COEFFICIENTS 

In  all  of  the  computations  with  the  quadratic  drag  formula  the 
following  sectional  viscous  drag  coefficients  were  used: 

Legs 

CD  =  0.50  for  legs  exposed  to  current 

CD  =  0.25  for  legs  which  are  directly  behind  another  leg 

Pontoons 

CD  =  0.45  for  pontoons  exposed  to  current 
CQ  =  0.45/2  for  pontoons  which  are  behind  other  pontoons 
CQ  =  0.20  for  the  smaller  10'  horizontal  member  between  the  pon¬ 
toons 


Vertical  Member 


Other  Small  Members 


Tendons  and  Risers 

CQ  =  0.60  at  free  surface  for  the  risers 

CQ  =  0.60  at  the  bottom  of  the  legs  for  the  tendons 

CD  =  1.20  at  1000-ft  depth  for  both  tendons  and  risers 

The  drag  coefficient  varies  linearly  from  top  to  1000-ft  depth 
and  is  constant  from  1000-ft  depth  and  down  to  the  bottom. 
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Figure  B-9.  Sway  displacement  for  maximum  regular  wave,  current, 
steady  and  unsteady  wind  and  wave  drift. 
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SWAY  RESPONSE:  (FEET)  SWRY  RESPONSE  (FEET  ) 

!0.0  -20.0  0.0  20.0  <0.0  -60.0  -10.0  -20.0  0.0  20.0  10.0 


SHAY  RESPONSE (FEET)  SWAY  RE5P0N5E  (  FEET  ) 


TIME l SECONDS) 

Run  B-6.  Sway  displacement  for  unsteady  wind  alone  (no  current, 
waves  or  steady  wind). 


TIME) SECONDS) 

Figure  B-7.  Sway  displacement  for  current  and  unsteady  wind  (no 
wave  or  steady  wind). 
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Figure  B-5. 


Sway  displacement  for  maximum  regular  wave, 
steady  wind  and  wave  drift  force. 


current, 
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Figure  B-2.  Sway  displacement  for  maximum  regular  wave  only 
(no  wind,  current  or  wa ve  drift  force). 
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Figure  B-3.  Sway  displacement  for  maximum  regular  wave  only, 

with  integration  up  to  disturbed  free  surface  level 
(no  wind,  current  or  wave  drift  force). 
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Table  B-2.  Surge  displacement  components 


Run  # 

Da 

xr  ft** 

x2,  ft*** 

1A 

14.4 

0.0 

0.0 

IB 

14.4 

0.0 

0.0 

2A 

0.0 

22.1 

0.0 

3A 

7.7 

22.1 

0.0 

4A 

16.1 

22.1 

0.0 

4B 

42.3 

22.1 

0.0 

5A 

70.3 

22.1 

0.0 

5B 

96.4 

22.1 

0.0 

6A 

0.0 

0.0 

40.9 

7A 

15.3 

0.0 

25.6 

7B 

15.5 

0.0 

warn 

8A 

0.0 

22.1 

15.0 

8B 

0.0 

22.1 

2.5 

9A 

70.5 

22.1 

WEM 

9B 

96.7 

22.1 

8.9 

*  XQ  =  steady  displacement 

**  Xj  =  first-order  (wave  frequency)  amplitude 

***  X2  =  second-order  (natural  frequency)  amplitude 

Wave  frequency  =  .450  rad/sec 
Natural  frequency  =  .600  rad/sec 
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smaller  for  the  interaction  between  the  low  frequency  motion  and  the  current 
than  for  the  high  frequency  motion  and  the  current  (see  Runs  4A  and  4B). 

Runs  8A  and  8B 

The  results  for  Runs  8A  and  8B,  with  the  maximum  regular  wave  and 
the  unsteady  wind,  are  shown  in  Figure  B-8.  These  runs  show  the  interaction 
between  the  low  and  high  frequency  motions.  Comparing  with  Run  6A,  it  can 
be  seen  that  the  low  frequency  motions  are  much  reduced  by  the  presence  of 
the  high  frequency  motions.  For  the  SAI  drag  formula,  the  change  in  the 
low  frequency  amplitude  is 

X*  =  15.0  -  40.9  =  -25.9  feet  (SAI  drag) 

and  for  the  quadratic  drag, 

X^  =  2.5  -  40.9  =  -38.4  feet  (quadratic  drag). 

Runs  9A  and  9B 

Figure  B-9  shows  the  results  for  Runs  9A  and  9B,  including  the 
maximum  design  wave,  steady  and  unsteady  wind,  current  and  wave  drift.  By 
comparing  with  Run  5B,  it  can  be  seen  that  the  low  frequency  motion  has 
little  effect  on  either  the  high  frequency  or  the  steady  surge. 


whereas  the  SAI  displacement  only  shows  a  small  increase.  It  is  believed 
that  the  quadratic  drag  fomrula  may  be  inadequate  in  predicting  this  inter¬ 
action  effect. 

Runs  5A  and  5B 

Figure  B-5  shows  the  results  for  Runs  5A  and  5B,  current,  steady 
wind  and  maximum  design  wave,  including  the  steady  second-order  drift  force. 

As  expected,  the  increase  in  steady  displacement  is  the  same  for 
both  drag  formulas, 

AXq  =  54.1  feet. 

The  increase  due  to  the  wave  drift  force  was  about  2.0  feet,  or  approximately 
9 %  of  the  primary  motion  amplitude. 

Run  6A 

Figure  B-6  shows  the  motions  for  the  unsteady  wind  only.  Run  6A. 

The  amplitude  of  the  motion  at  the  surge  natural  frequency  is 

X2  =  40.9  feet. 

Runs  7A  and  7B 

Figure  B-7  gives  the  results  for  Runs  7A  and  7B  with  current  and 
unsteady  wind.  Comparison  with  Run  6A  shows  the  change  in  the  low  frequency 
motion  due  to  the  current.  For  the  SAI  drag  formula,  this  change  is 

AX^  =  25.6  -  40.9  =  -15.3  feet  (SAI  drag) 

and  for  the  quadratic  drag  the  change  is 

AXC2  =  20.7  -  40.9  =  -20.2  feet  (quadratic  drag). 

It  can  be  seen  that  the  discrepancy  between  the  two  methods  is  considerably 


